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Simuiationsvorrichtung (1), die die Momentenmethode 
verwandet und mit der die intensitat eines elektromagneti- 
schen Feldes etc. im Zehbereich mit hoher Geschwindigkert 
simuliert werdan kann. Die Simuiationsvorrichtung (1) um- 
faSt TransformationsmitteJ (11) rum Auffinden eines Fre- 
quenxspektnims, indem eine zeitliche Datanfolge einer 
Welle nquelle untertoitt und darauf die Fourier-Transforma- 
tion angewendet wird, Berechnungsmittel (12) zum Berech- 
nen der gegenseKigen Impedanz fur eine Abtastfrequenz. 
Ermittaln einer Naherungsgleiehung fur die gegenseltige 
Impedanz (Z«) anhand der berechneten gegensertigan Impe- 
danz und der Abtastfrequenz, und zum Berechnen der 
gegenseitigen Impedanz fur bestimmte von den Transforma- 
tionsmitteln (11) vorgegebene Frequenzen anhand der Nahe- 
rungsgleichung, Simulationsmittel (13) zum Ermitteln des in 
jedem Element flieSenden Stromspektrums durch Anwen- 
dung der Momentenmethode anhand der durch die Berech- 
nungsmittel (12) berechneten gegenseitigen Impedanz und 
des von den Transformationsmitteln (11) vorgegebenen 
Frequenzspektrums, und zum Ermrttem eines Spannungs- 
spektrums, eines elektrischen Feldspektrums und eines 
magnetischen Feldspektrums anhand des Stromspektrums, 
und Inverse-Transformationsmrttel (U) zum Anwenden der 
inverse n Fourier-Transformation auf die von den Simula- 
tionsmitteln (13) ermittelten Spektren und zum Ausgeben 
des Ergebnisses. 
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Beschreibung 

Die voriiegende Erfmdimg betrifft eine Simn la tionsvorrichtung and an Simulattonsverfahren Simulieren 
der Intensitat ernes von einer elektronischen Vorrichtung abgestrahlten elektromagnetischen Feldes etc mit 
5 Hilfe der Momentenmethode, insbesondere betrifft die vorliegende Erfmdimg eine Simulationsvomchtung und 
em die Momentenmethode anwendendes Simula tionsverfahren, womit die Intensitat ernes elektromagnetischen 
Feldes etc im Zeitbereicfa mh einer hohen Gesc&windigkett simuliert werden fcarm 

Die Abstrahlong ongewunschter RadioweHen oder von Rauschen durch elektronische Vorrichtungen wurde 
auf bestimmte zulassige Pegel beschrankt. Diesbezugiich wurden von den emzelnen Staaten strenge Normen 

10 oder Vorschriften auf gestellt 

Znr Erfullung derartiger Radio wellenvorscbrift en werden verschiedene Gegenmessungen, wie z. B. bei der 
Ahschirmungs- oder FBtertechnik, angewendet. Werden derartige Gegenmessungen angewendet, ist die Ent- 
wicklung einer Simulationstechnik erforderiich, um quantitativ den Verringerungsgrad der Radiowellen oder 
des Rauschens simulieren zu kdnnen. 

15 Ebenso wurde die Beeinflussung elektronischer Vorrichtungen durch von anderen elektronischen Vorrichtun- 
gen abgestrahlte Radiowellen oder Rauschen auf einen bestimmten zulassigen Pegel beschrankt, der hdher als 
der oben genannte festgelegte Pegel ist Auch diesbezugiich wurden von den emzelnen Staaten strenge Normen 
oder Vorschriften auf gestellt. 

Um diese RadioweUenvorschriften erfuDen zu konnen, ist die Entwicklung einer Simulationstechnik erforder- 

20 ftch, so dafi die Grunde dafur ermittelt werden konnen, warum elektronische Vorrichtungen ungewunschte 
Radiowellen oder ungewunschtes Rauschen abstrahlen und warum Fehlfunktionen von elektronischen Vorrich- 
tungen aufgrund der Radiowellen oder des Rauschens auftreten. 

Zur Entwicklung der zuletzt genannten Simulationstechnik ist eine Vorrichtung zur Berechnung der intensitat 
eines elektromagnetischen Feldes erforderiich, womit die Intensitat eines sich -rritiirh andemden elektromagne- 

25 tischen Feldes simulieren werden kann. Eine derartige Vorrichtung zur Berechnung der Intensitat eines elektri- 
schen Feldes, um die Intensitat eines sich zertfich andemden elektromagnetischen Feldes berechnen zu konnen, 
ist jedoch bisher noch nicht in die Praxis umgesetzt worden. Nachfolgend werden die Grunde hierfur erlautert. 

Die Intensitat eines von einem betiebig gefbrmten Objekt abgestrahlten elektromagnetischen Feldes kann auf 
einfache Weise mit Hilfe bekannter theoretischer Gleichungen berechnet werden, falls der in jedem Teil des 

30 Objekts flieBende elektrische und magnetische Strom bekannt sind. Die Werte des elektrischen und magne- 
tischen Stromes konnen bei vorgegebenen Rgnrfh^m^myn durch Losen der elektromagnetischen Max- 
wefl'schen Gleichungen erhalten werden. 

Ein Verfahren zum Losen dieser Gleichungen ist die Momentenmethode. Die Momentenmethode ist eines der 
Verfahren zum Losen von aus den elektromagnetischen Maxwell'schen Gleichungen abgeleiteten Integrations - 

35 gleichungen, wo bei die Momentenmethode der elektrischen und rnagnetischen Strom berechnet, indem das 
Objekt in kleine Elemente aufgeteilt wird, so da6 mit der Momentenmethode jedes beEebige dretdimensional 
geformte Objekt behandeh werden kann. Eine die Momentenmethode beschreibende Dnickschrift ist beispiels- 
weise "Sinusoidal Reaction Formulation For Radiation And Scattering From Conducting Surface* H. N. Wang, J. 
H. Richmond und M. C Gureath, IEEE TRANSACTIONS ANTENNAS PROPAGATION, voL AP-23, 1975. 

40 Andererseks ist eine Analyse im Zehbereich erforderiich, um die Grunde dafur ermitteln zu tannin, warum 
eine elektronische Vorrichtung ungewunschte Radiowellen oder ein ungewunschtes Rauschen abstrahlt und 
warum bei einer elektronischen Vorrichtung eine Fehlfunktion aufgrund von Radiowellen oder Rauschen 
auftritt, da die meisten Fehlfunktionen von elektronischen Vorrichtungen durch pulsahnHches Rauschen verur- 
sacht werden. Des weiteren werden haufig Fehlfunktionen von elektronischen Vorrichtungen durch einen 

45 abnormalen Betrieb von in den Vorrichtungen enthattenen Bauteilen, wie z. R ICs, verursacht. Diese Bauteile 
mussen zeitfich uberwacht werden, um einen abnormalen Betrieb dieser Bauteile feststeflen zu konnen. 

Zwei nutzliche Verfahren zur Analyse im Zehbereich sind die finite Elementenmethode und die finite Diffe- 
renzenmethode. Obwohl mit der finken Elementenmethode und der fmiten Diff erenzenmethode eine Analyse im 
Zehbereich moglich ist, konnen mh diesen Verfahren elektronische Vorrichtungen, die eine Viehahl von Struk- 

50 turen wie z. R Obertragungslertungen, Kabel und Gehause usw. aufweisen, nur schwer behandelt werden. 

Der Grund hierfur ist die Tatsache, dafi bei der fmiten Elementenmethode und der fmiten Diff erenzenmetho- 
de das Objekt und der das Objekt umgebende dreidimensionale Raum unterteih werden mussen. Werden kleine 
Teiie, wie z. R Kabelendbereich, sehr fein unterteih; wird die Anzahl der Untertetlungen enorm hoch, da der das 
Gehause und die Kabel umgebende Raum sehr groB ist, was jedoch zu einer Oberbelastung der internen 

55 Spetcherkapazxtaten des Computers fuhrt Wird hingegen fur die Kabel, das Gehause oder fur andere Strukturen 
lediglich eine grobe Unterteilung durchgefuhrt, kann der EinfluS von Kabelendbereichen unmogiich analysiert 
werden, die jedoch in dem Gesamtmechanismus eine bedeutende Rolle spielen und wichtige QueOen fur 
ungewunschte Radiowellen oder ungewunschtes Rauschen darstellen. 

Des weiteren wird als Koordinatensystem fur die finite Elementenmethode und die finite Diff erenzenmethode 

60 gew6hnlich das kartesische Koordinatensystem verwendeL Die Kabel und Kabelendbereich e, die in dem Ge- 
samtmechanismus eine bedeutende Roile spielen, umfassen jedoch zylindrische Elemente, wahrend das Gehause 
der Vorrichtung betiebig geformt sein kann. Somit entsteht das schwer zu Idsende Problem, wie das zu analysie- 
rende Objekt unterteilt werden solL 

Diese Probleme bestehen bei der Momentenmethode nicht, und die Momentenmethode ist zur Behaudlung 

65 von elektronischen Vorrichtungen, die eine Vielzahl von Strukturen wie z. B. tJbertragungsleitungen, Kabel oder 
Gehause umfassen, best ens geeigneL 

Der Grund hierfur ist die Tatsache, dafi die Momentenmethode eine Art Begrenzungselementenmethode ist, 
bei der lediglich die Begrenzungsoberflache zweidimensional unterteilt werden tnuB. Des weiteren kann der 
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UnterteHungsgrad ^^1^^^^^ ^emt/^odt^dt^S 

SSX^weS^SXS^ -es vertraghchen Zci.ahn.cns 

wW erfindangsgemae dnrd. cine Sim^atK,^r^tung nach den Anspruchen 1, 2 oder 10 
undein Simulationsverfahren nach den Anspruchen 12, 13 oder « gelost. _ , 

""fE^I^^Spt^p und den Aufbau der voriiegenden Erfindun& _ _ , m „ 

£!; j Wt erne St Darsteuung des Prinzips and des AufW der vorbegenden Er&idung, 
* ^SS ^VeDaiSellunE zur ErUuteruiig der Ermittiung einer Nihenmgsgleidnm& 

r^SA^sSn^^^^^^ 

Fte 6 zeiet eme vierte DarsteUung zur Erlautenmg der Ermitthmg der Naheruiigsg leichung, 
£!■ 7 Darstelh.nl zur Erlautenmg der Ennitthing der NahenmgsgUaetang. 

K* 8 " S^Se D^XngzurEriiuterungder Ermitthmg der ^cnir^glcrfiang, 
Fte 9 zeto etae siebte DarsteUung zur Erlautenmg der Ennitthingder NSheruiigsgleictang. 

m zeto^e ach«e D^elhirfg zur Erlautenmg der Ermittiung der Naherungsglcchung, 
S nrilB^d nSSne^e DarsteUungen zur ErUutening der ^^^ a ^ e,da0g - 
irf" 15Teirtefae zehnte DarsteUung zur Erlimerung der Ermittiung der Naherungsglcchung, 
S USdSSSSta DarSeuungen zur Eriauterung der ErmMung der Nlh*«.tfacfa»g. 
Fte 4 zeta eine zw61fte DarsteUung zur Erlautenmg der Ermittiung der Naherur^glcKiung 
Fe 15 zcta etae dreizebnte DarsteUung zur Erlautenmg der Ermittiung der Naherungsg etchung, 
Si' ^"Se^e^etaiteDarstellunizurErUutenmgderErm 

Fta 17 "Setae Bnfeehnte DarsteUung ror Erlautenmg der E^^ d ^ Naherung S gl«chumj 
Fg'. S cine sedLhnte DarsteUung zur Erlauterung der Ermittiung der Naherungsgle.chung. 
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Fig. 1 9 zeigt em Ausfunrungsbeispiel der voriiegenden Erfmdung, 
Fig. 20 zeigt ein Ausfuhrungsbeispiel fur das Kngeben einer Wellenquelle, 
Fig. 21 zeigt eine DarsteUung zur Eriauterung des Eingebens einer WeilenqueUe, 
Fig. 22A und 22B zeigen Darstelfamgen einer Wellenquellen fur etnen Impuistest, 
Fig. 23A und 23B zeigen weitere Darstellungen der Wellenquelle fur den Impubtest, 
Fig. 24 zeigt ein Ausfuhrungsbeispiel zur Berechnung ernes Naherungswertes, 
Fig. 25 zeigt ein Ausfuhrungsbeispiel zur Berechnung eines OberweUenstromes, 
Fig. 26 zeigt Simulationsgleichungen gemaB der Momentenmethode, 
Fig. 27 zeigt ein Ausfuhrungsbeispiel zur Berechnung eines Stromes, 
Fig. 28 zeigt ein Aiisfuhnmgsbeispid zur Berechnung einer Spannung, 
Fig. 29 zeigt ein Ausfuhrungsbeispiel zur Berechnung eines elektromagnetischen Feldes, 
Fig. 30 zeigt eine erste Darstellung zur Eriauterung der Ermittlung von einer zwischen Leitern anliegenden 
Spannung, 

- Fig. 31 A, 3 1 B und 31 C zeigen zweite Darstellungen zur Eriauterung der Ennhthmg der zwischen den Leitern 
anliegenden Spannung, 

Fig. 32 zeigt eine erste Darstellung zur Eriauterung der Berechnung der Kopphmg einer Obertragungsleitung, 
Fig. 33 zeigt eine zweite Darstellung zur Eriauterung der Berechnung der Kopplung der Obertragungsleitung, 
Fig. 34 zeigt eine drrtte Darstellung zur Eriauterung der Berechnung der Kopplung der Obertragungsleitung, 
Fig. 35 zeigt eine vierte DarsteUung zur Eriauterung der Berechnung der Kopplung der Obertragungsleitung, 
Fig; 36A und 35B zeigen funfte Darstellungen zur Eriauterung der Berechnung der Kopplung der Obertra- 
gungsleitung, 

Fig. 37 zeigt eine sechste Darstellung zur Eriauterung der Berechnung der Kopplung der Obertragungslei- 
tung, 

Fig. 38 zeigt eine siebte DarsteUung zur Eriauterung der Berechnung der Kopplung der Obertragungsleitung, 

Fig. 39 A und 39B zeigen erste Darstellungen zur Eriauterung der Simulation, 

Fig. 40 zeigt eine zweite Darstellung zur Eriauterung der Simulation, 

Fig. 41 zeigt Simulationsgleichungen gemaB der Momentenmethode, 

Fig. 42 zeigt eine erste Darstellung zur Eriauterung einer Strahlungsanalyse einer Antenne, 

Fig. 43 zeigt eine zweite Darstellung zur Eriauterung der Strahlungsanalyse der Antenne, 

Fig. 44 A und 44B zeigen drhte Darstellungen zur Eriauterung der Strahlungsanalyse der Antenne, und 

Fig. 45 zeigt eine Darstellung zur Eriauterung der Momentenmethode. 

Bei der Berechnung der Strahlungsintensitat eines elektromagnetischen Feldes mit HDfe der Momentenme- 
thode wtrd, wie durch das FluBdiagramm von Fig. 45 dargestellt, beun Lesen von Konfigurationsmformationen 
der zu simnlierenden gitterartjgen elektronischen Vorrichtung eine bisher noch nicht verarbehete Frequenz aus 
den auszuwertenden Frequenzen ausgewahlt, und fur diese Frequenz werden nut Hilfe eines bestimmten 
Berechnungsverfehrens die gegensehige Impesrlanz, die gegenseitige Adnuttanz und die Wechselwirkung zwi- 
schen den grtterartigen Elementen ermrtteh. Die somit ermhtelte gegenseitige Impedanz usw. und die durch die 
Kjonngurationsinfbrmauonen spezuxrierte Wellenquelle werden in den Simulationsgleichungen gemaB der Mo- 
mentenmethode eingesetzt, und durch Ldsen der Simulationsgleichungen wird der in jedem Element flieBende 
elektrische und magnetische Strom ermittelt Diese Strdme werden anschliefiead zur Berechnung der Strah- 
tungsintensrtat an jedem Untersuchungspunkt eines elektromagnetischen Feldes verwendet. Die Berechnung 
jeder gegenseitigen Impedanz, gegenseitigen Adnuttanz und jeder Wechselwirkung benotigt jedoch eine lange 
Zeh, und diese Werte mussen zudem fur jede einzelne Frequenz des Frequenzbereiches ermittelt werden. Wie 
beretts erwahnt, kann somit die Berechnung in der Praxis nicht innerhalb einer kurzen Zeuspanne durchgefuhrt 
werden. 

Insbesondere betragt zwar, wie bereks beschrieben, die fur die Ldsung der Simulationsgle tchungen gemaB der 
Momentenmethode bendtigte Zeit fur eine einzige Frequenz einige Minuten; die Berechnung der gegenseitigen 
impedanz, der gegenseitigen Admittanz und der Wechselwirkung benotigt jedoch einige Stunden. Da diese 
Werte fur jede Frequenz des Frequenzbereiches berechnet werden mussen, kann diese Berechnung nicht 
innerhalb einer praktisch vertraglichen kurzen Zeitspanne durchgefuhrt werden. 

Die voriiegende Erfmdung ermoglicht hingegen eine sehr schnelle Simulation der Intensitat eines sich zeitCch 
andernden elektromagnetischen Feldes usw. Nachfolgend wird die Erfmdung detailliert beschrieben. Es sei 
darauf hingewiesen, daB bei der erfmdungsgemaoen Berechnung der Strahlungsintensilat eines elektromagneti- 
schen Feldes mit Hilfe der Momentenmethode vorzugsweise erne Technikfur die Hochgeschwmdigkeitsrealisie-- 
rung der Momentenmethode zum Einsatz kommt, die bereits zuvor von den Erfrndern der voriiegenden Erfm- 
dung erhmden und beschrieben worden ist (vgL Japanische Patentamneldung 7-298062). Es wird daher zunachst 
die von den Erfrndern der voriiegenden Erfmdung erfundene Hochgeschwmdigkettsrealisierung der Momenten- 
methode erlautert 

Hinsichtlich der Hochgeschwmdigkeitsrealiaerung der Momentenmethode wird nachfolgend die Frequenz 
mit f, die Wellenzahl mit k (=2m7c, wobei c die Uchtgeschwindigkeit angibt) und der grundsatzliche Abstand 
zwischen einzelnen Elementen mit r 0 bezeichnet 

Die zwischen einem Element i und einem Element j, die nach der Momentenmethode aufgeteth worden sind, 
auftretende gegenseitige Impendanz Zij kann durch folgendes Frequenz-Porynom riaherungsweise ausgedruckt 
werden: 

Zij - e-i^Rao + a t P + a 2 f* + ajf 5 + a4f» + . . .) + jCbof 1 + b,f + fcf 3 + baf 5 + b^f 7 + ...)J 

Die gegenseitige Admittanz Yij zwischen den Elementen i und j kann ebenfalls durch folgendes Frequenz- 
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Polynom naherungsweise ausgedriickt werden: 

Yy « e-i^Rao + a,P + a 2 f* 4- ajf 6 + tuP + ."..) + jW 1 + W + W 1 + W* + W + »-)]. 

Die Wechseiwirkung zwischen den Elementen i und j kann schlieBlich durch folgendes Frequenz-Poiynom 
riaherungsweise dargesteUt werden: 

Bij « e-i^Rao + a t P + a 2 f* + ajf 5 + suf 8 + ...) + j(bof + W 3 + Iwf 5 + hat 7 . + M 9 » -ft 

Zunachst werden die gegenseitige Impedanz Zi h die gegenseitige Admittanz Yy und die Wecfaseiwirkung fig 
mit Hilfe eines genauen Berechnungsverfahrcns fur unterschiedliche Abtastfrequenzen berechnet Die sicn 
daraus ergcbenden Berecfanungswerte werden in die obigen Nahenmgsglriduingen eingesetzt, urn auf diese 
Weise die Simulationsgleichungen znm Ermitteln der Koeffizienten a p und bp fur jeden werteren Berecnnungs- 
wert aufzustellen, und diese Simulationsgleichtmgen werden nachfolgend gelost, urn somit die Koeffiaenten ap 
und b, zu ermitteln. Wird anschiieBend eine von den zunachst verwendeten Abtastfrequenzen abweicnende 
Frequenz vorgegeben, werden die gegenseitige Impedanz Z,> die gegenseitige Admittanz Yg und die Wecasel- 
wirkung Bij jeweils mit Hilfe der obigen Naherungsgieichungen unter Verwendung der zuvor ermrtteken 
Koeffizienten a p und bp berechnet Auf diese Weise kann die Momentenmethode mit hoher Geschwmdigkeit 
realisiert und implementiert werden. , t _ . . . . , 

Essei darauf hingewiesen, daB die Anzahl der zunachst verwendeten Frequenzen davon abhangig^t, wieviele 
Koeffizienten ap und bp der Naherungsgieichungen ermittelt werden sollen, Werden beispielsweise funf Abtast- 
frequenzen verwendet, konnen die Koeffizienten ao bis 34 und bo bis b4 ennittek werden. 

Wie nachfolgend eriautert wird, wird erfindungsgemSB ein Verfahrenvorgesddagen beidem eine sichzeithch 
andernde WeUenquefle zunachst in den Frequenzbereich transformiert, anschiieBend die Momentenmethode 
zur Simulation des Strom- oder Spannungsveriaufes der einzemen Hemente sowie ^f^^^^f 1 ^ 
elektromagnetischen Feldes in dem Transfonrianons-Frequenzbereich verwendet und scMeBfich die Sunuia- 
tionswertein den Zehbereich wieder zoriick- bzw. invertransf onniert werden. Wird dabei die oben beschnebene 
Hochgeschwindigkeitsrealisiening der Momentenmethode eingesetzt, kann die gegenseitige Impedanz die 
gegensertige Admittanz Y i3 und die Wechseiwirkung Bij sehr schneD berechnet werden, so daB die Momenten- 
methode mit hoher Gescfawindigkeh ansgefuhrt werden kann. . . 

Auf diese Weise kann innerhaib kurzer Zeit festgesteUt werden, warum eme elefaroniscne Vornchtung 
ungewfinschte RadioweUen oder ungewunschtes Rauschen abstrahlt und warum eme Fehlfunktion einer elek- 
tronischen Vornchtung aufgrund von RadioweUen oder Rauschen auftritt. 

Das Prinzip sowie der grundsatzfiche Aufbau der vorliegenden Erfindung werden nachfolgend unter Bezug- 
nahm e auf Fist 1 und 2 erl&utert. . . . 

In diesen Figuren bezeichnet das Bezugszeicfaen 1 eine erfmdungsgemaBe Simuiationsvomchtung. Diese 
Vornchtung teilt eine elektronische Vornchtung in emzelne Bemente auf, berechnet fur erne vorgegebene 
Frequenz die zwischen den Hemeaten auftretende gegenseitige Impedanz, gegemer&ge Admittanz und Wecn- 
selwUung und simuUert gemaB der Momentenmethode aufgrund der berechneten Werte die m jedem Element 
ffieBenden elektrischen und magnetischen Strome sowie die den Hementen entsprechende Weflenqueue. 

Die in Fig. 1 dargesteflte erfmdungsgemaBe SimiiUtwnsvorridrtung 1 umfafit Verwaltmigsmittel 10, Transfor- 
mationsmittelll,Berechnimgsmitteil2,Suniilad . , 

Die Verwaltungsmittel 10 verwaken die Krafiguralionsmf ormationen der zu analysierendeneletomschen 
Vornchtung, welche zur Simulation in eine Vielzahl einzemer Hemente aufgeteflt worden 1st. Die Transforma- 
tionsmittell 1 teUen die zehabhangigen Daten der WeflenqueUe in eine zeMiche Datenf olge auf tmd ffihren die 
Fouriertransfonnation durch, urn auf diese Weise das entsprechende Frequenzspektrum zu erhalten. Die Trans- 
fonnationsmittel tl legen das gesamte Frequenzspektrum oder ein Teil des Frequenzspektrums bzw. die 
entsprecfaenden Frequenzen als Simulationsobjekt fest tv«„^™ 

Die Berechnungsmhtel 12 umfassen Auswahlmittel 15, urn unter Berucksichtigung der von den Transforma- 
tionsmittem 11 festgdegten Frequenzen eine Abtastfrequenz auszuwahlen, sowie erste Ben^ungsimttel 16 
zur Berechnung der gegensehigen Impedanz, der gegensehigen Admittanz und der Wechseiwirkung fur che von 
den Auswahlmittem 15 ausgewafahe Abtastfrequenz. Des werteren weisen die Berechnungsmmel 12 Erzeu- 
Eunesinittel 17 auf, urn eine Naheningsgteichung fur die gegenseitige Impedanz auf Grundlage der Abtastfre- 
quenz und der von den ersten Berechnungsmittem 16 fur die Abtastfrequenz berecfaneten gegenseitigen Impe- 
danz, eine Naherungsgleichung fur die gegenseitige Admittanz auf Grundlage der Abtastfrequenz und der von 
deTersten Berechniingsmitteln 16 fur die Abtastfrequenz berechneten gegenseitigen Admittanz und eme 
Naherungsgleichung fur die Wechslwirkung auf Grundlage der Abtastfrequenz und der von den ersten Berech- 
nunesmitteui 16 fur cue Abtastfrequenz berechneten Wechseiwirkung zu erzeugen. Die Berechnungsmittel 12 
umfassen schlieBlicfa auchzweite Berechnungsmittel 18, urn fur weitere von den Transformationsmitteln vorge- 
gebene Frequenzen mit Hilfe der durch die Erzeugungsmittel 17 erzeugten Nahenmgsgleichungen die gegensei- 
tige mipedanz,cfiegegemeitigeAdnuttanz und dieWechseiwiriorngzuberechnen. ^ 

Die Shnu^nsmittel 13 ermhtem aufgrund der von den Berechnungsmittem 12 berechneten Werte fur die 
gegenseitige Impedanz, gegenseitige Admittanz und Wechseiwirkung sowie aufgrund des von den Transforma- 
lonsm\ttem 11 vorgegeblnen Frequenzspektrums mit Hilfe der Momentemriethode das m den canto 
Elen^nten flieBenden elektrische und magnetische Stmmspektrum. Gleichzeitig ermitteln d« |SimuUnonsmittel 
13, falls eewunscht abhangig von dem elektrischen Stromspektrum und magnetischen Stromspektrum das 
entsprechende Spannungs-^ elektrische Feld- und magnetische Feldspektrum und legen fest welche gewunsc^ 
ten Fre^enzen des ernStelten Spektnims ausgegeben werden sollen. Die Inverse-Transformationsmrttel 14 
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transformieren mh Hilfe der inversen Fouriertransformation das von den Simulationsmitteln 13 ausgegebene 
Frequenzspektrum in den Zeitbereich und geben das Ergebnis axis. 

Die Simulatkrasmittel 13 fuhren somit gegebenenfalls die Simulation durch, urn abhangig von einem entspre- 
chenden Befehl einer Bedienperson ledigEch das elektrische und magnetische Stromspektrum zu ermhteln. In 
5 dies em Fall konnen mit Hilfe der in Fig. 1 gezeigten J/M-Mhtel 19 die entsprechende Spanming and das 
entsprechende elektrische and magnetische Feld im Zeitbereich abhangig von dem im Zeitbereich vorliegenden 
und von den In verse-Transf ormatio nsmitteln 14 ausgebenen elektrischen und magnetischen Strom ermitteh 
werden. 

Nachdem die Transfonnationsmittel 11 dorch Unterteilen der in zerdicher Folge vorliegenden Daten der 
io WeQenquelle und Durchfuhren der Fouriertransformation das Frequenzspektrum erhatten und das gesamte 
Frequenzspektrum oder einen Teil davon sowie die entsprechenden Frequenzen fur die nachfolgende Verarbei- 
tung festgelegt haben, wahlen die Auswahlmittel 15 der erfmdungsgern a Ben Simulationsvorrichtung 1 unter 
Berucksichtigung der von den Transformationsmitteln 11 vorgegebenen Frequenzen eine Abtastfrequenz aus, 
und die ersten Berechnungsmittel 16 berechnen mit Hilfe ernes genauen Aigorithmus die gegenseitige Impedanz, 
15 die gegenseitige Admittanz and die Wechserwirkung fur die von den Aaswahlmittem 15 ausgewahlte Abtastfre- 
quenz. 

Anschliefiend ermittem die Erzeugtmgsmittel 17 die Simulationsgleichungen, wobei die Koeffizienten ap and 
bp unbekannte GroBen darstellen, indem die gegenseitige Impedanz Zij beispielsweise durch folgenden Aus- 
druck naherungsweise beschrieben wird 

Zij - e-^Rao + aif 2 4- a 2 f* + asf 5 + suf* + ...) + j'Cbof 1 + bif + bsf 3 + bsf 5 + b*! 7 + ...)] 

(wobei f die Frequenz, k die Wellenzahl and ro den Abstand zwischen den FJementen i und j darsteUt) und die 
Abtastfrequenz und die von den ersten Berechnungsmitteln 16 berechnete gegenseitige Impedanz in diese 
25 Nanerungsgleichong eingesetzt und die sich daraus ergebenden Gfeichungen gelost werden. 

Des weheren ermhteln die Erzeagungsmittel 17 die Simulationsgleichungen, wobei die Koeffizienten a p and 
bp unbekannte GroBen darstellen, indem die gegenseitige Admittanz Yij beispielsweise durch folgenden Aus- 
druck naherungsweise beschrieben wird 

30 Yy = e-^Rao + aif 2 + a 2 f* + ajf 6 + *4* + ...) + jfbof 1 + b t f + b^f 3 + bjf 5 + b^f 7 + ...)] 

(wobei f die Frequenz, k die Wellenzahl und ro den Abstand zwischen den FJementen i und j darsteUt) und die 
Abtastfrequenz und die von den ersten Berechnungsmitteln 16 berechnete gegenseitige Impedanz in diese 
N3herungsgleichang eingesetzt und die sich daraus ergebenden deichnngen gelost werden. 
35 Schfiefilich ermhteln die Erzeugungsmhtel 17 auch die Simulationsgleichungen mit den Koeffizienten ap and 
bp als Unbekannte, indem die Wechselwirkung Bij beispielsweise dorch folgenden Ausdruck naherungsweise 
beschrieben wird 

Bij - e-^Rao + a^ + a 2 f* + asf 5 + a*? + . . .) 4- j(bof + biF + b 2 f 5 + b 3 f 7 + b4f 9 + -..)] 

40 

(wobei f die Frequenz, k die Wellenzahl und r 0 den Abstand zwischen den FJementen i und j darsteUt) und die 
Abtastfrequenz and die von den ersten Berechnungsmitteln 16 berechnete Wechselwirkung in diese Nanerungs- 
gleichung eingesetzt und die sich daraus ergebenden Gieichungen gelost werden. 

Nachdem auf diese Weise die Naherungsgleichungen fur die gegenseitige Impedanz, die gegenseitige Admit- 

45 tanz und die Wechselwirkung bestimmt worden sind, berechnen die zweiten Berechnungsmittel 18 fur jede 
weitere von den Transformationsmitteln 11 vorgegebene Frequenz die entsprechende gegenseitige Impedanz, 
gegenseitige Admittanz und Wechselwirkung durch Anwendung dieser drei Naherungsgleichungen (fiber den in 
Fig. 1 dargestellten Pfad LI). Nach Berechnung dieser Werte ermitteln die Simulationsnuttel 13 mit HOfe der 
Momentenmethode aufgnmd der berechneten gegenseitigen Impedanz, gegenseitigen Admittanz und Wechsel- 

50 wirkung sowie auf grand des von den Transformationsmitteln 1 1 f estgelegten Frequenzspektnims das Spektrum 
des in jedem Element flicBenden elektrischen und magnetischen Stromes (fiber den in Fig. 1 dargestellten Pfad 
L2). 

Anschliefiend bestimmen bzw. spezifizieren die Simulationsmittel 13 dasjenige Spektrum des elektrischen und 
magnetischen Stromes der ermittehen Stromspektren, fur welches die Simulation gewunscht wird bzw. welches 

55 fur die weitere Verarbeitung verwendet werden sott, sowie das Spektrum der Spanming und des elektrischen 
sowie magnetischen Feldes, fur welches die Simulation gewunscht wird, und geben diese Spektren aus bzw. 
bestimmen diese Spektren fur die Ausgabe. Die Inverse-Transfonnanonsmittel 14 transformieren die ausgege- 
benen Spektren mh Hilfe der inversen Fouriertransformation in den Zeitbereich und geben das Ergebnis aus. 
Auf diese Weise werden durch die in Fig. 1 gezeigte Simulationsvorrichtung der vorliegenden Erfindung die 

go gegenseitige Impedanz, gegenseitige Admittanz und Wechselwirkung mh hoher Geschwindigkeit durch Anwen- 
dung der Naherungsgleichungen berechnet. Somit kann bei einer sich zeitlich andemden Wellenqueflc diese in 
den Frequenzbereich transformiert und die Intensitat des elektro magnetischen Feldes etc. mit HOfe der Momen- 
tenmethode in dem Frequenzbereich simcliert werden, wobei anschHeSend die Simulationswerte in den Zeitbe- 
reich zurucktransf ormiert werden. 

65 Mh Hilfe der erfindungsgemaBen Simulationsvorrichtung 1 kann daher der EinfluB einer Welienquelle auf 
eine elektro nische Vorrichtung im Zeitbereich simuliert werden, so dafi analysiert werden kann, warum eine 
elektronische Vorrichtung ungewunschte Radiowellen oder ungewunschtes Rauschen abstrahlt und warum in 
einer elektronischen Vorrichtung eine Fehlfunktion aufgrund von Radiowellen oder Rauschen auftritt. 
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Die in Kg. 2 gezeigte Simutauonsvorrichtong 1 der vorfegeaden E^g^^ff^ Verwahungs- 

M oauiiiuu«*»6 « ;»4»fTi F1«nent fliefienden elektrischen and maguetiscfaen Strom. Die *requenz 

S?£SSta£ ^m'SSSSSS Special CoZhtee on Radio Interference) 

^Ber^S^m'iSnen den elektrischen und magnetiscnen Strom im Zeitbcrei^ Die Ermjtt- 
lun^ttdMe^rXarfgr^ddes yon den BerechnungsnAtein 23 berectoetea und un Zeitbereicfc vorhe- 
Se^ek^Tn ^d Sdschen Stromes den gewOnscbten Strom, die gewunschte Spannnng bzw. das 

^W^ulteTSS So^^Weaenquene verwendet, teilen die Trennmittel 2. 
/wS5ffi% ^SgCeritoAmgsgemaBen Stoulationsvorrichtung 1 mdrei WeHenquel- 
J^JSweus anter^eK'KdTe TrSgersigaalfrequenz und Modnktioi^Enaifreque^ defimerte 
FTeo^'i^V a^TolrsStionSttel 22 S Sn durch Anwendung der Mcmemenmethode auf die 
dnrchSeSeTM aufgeteihen WeUenqueUen den in jedem Element fheBenden elektnschen und magne- 

^SSdtllnen auf Gnmdlage der Simulationsergebnisse den dem durch die Simuk- 
riJ^ftte^ta^qoenSerlS^mittelten elektrischen bzw. magnetischen Strom entsprechenden elefcn- 
uonsmtftel 22 un ZT^Tn MAereicb, AnschlieBend enmttein die Ermitdungsmittel 23 durch gemem- 

de^MAeSyorfegende elektrische und magne&che Strom vor deren Kombmation ausgewertet und erst 

^ttSS^ ^ modufiertes Trigersigna! ansgebende WeUenqaeUe verwendet^ wurde bei der 
in^l^fe^sS^^onsvomchttmg. bei der das Frequenzspektrum durch UnterteUung derm zerthcher 
?JE?«jfe^«]« Mote WeUenqueUe und anschlieBender Anweiidur«derFounertraiisfo™nongewon- 
FoIgevorBegenden Dates der J^™*™*^^ hoch wer den. Daher wird das im Zusammenhang nut Fig. 1 
bt^ebete g^y£2foteS£^£*" WeQenqueUe yerwende,, die em modules Tragers> 
Zta«iMkta FaU erfindungsgemaB die in Kg. 2 gezeigte Snnulator^ornchnrng 
^en^So^ sSrTjeSem Element ffieBende elektrische und magnetiscte Strom auf emfache Wese 

^SS^^^^^^r unter Bezugnahme auf verschiedene Ausrahrungsbei- 

'^Si wird ein Verfahren vorgescblager, mit dem -™^^>^^™^J^ 
un^An^Sfder Momentenmethode die Intensity eines elektromagnetischen Feldes auBerst schnell 

^aSluMM^ wird die zwischen einem Element i und einem Hement j auftretende gegenserdge 
ImpedaraZijwiefolgtr^enmgsweisebeschrieben: 

Die zwischen den Hementeni und j auftretende gegenseMgeAdrm^ 
naherungsweise beschrieben: 

Yy . e-fOR* + a,P + a.f* + aaP + a,? + ...) +KboP + b,f + b,P + b3p + brf 7 + -» 

Pie Wechsetwirkung By (elelctrischer Stromdipol gegenuber magnetischer Stromdipol) zwischen den Hemen- 
ten i und j wird folgendermaBen naherungsweise beschneben: 

Bjj _ e -i«[(ao + a,P + a 2 f« + ajf 6 + 84? + ...) +KW + b,P + b 2 P + W + M» + ...)} 

DabeibezeicbnetfdieFrequer^kdieWeUeroaU^^ 

den grundsatzlichen Abstand zwischen den Elementen l und j. . . • 7 _ _ it mfl , ein „ 

Nacnfolgend wird obige Naherungsgleichung zur Berechnung der gegensemgen Impedanz Z, nut Hflf e ernes 

^^n^^^mssElc^ fur die gegenseitige Impedanz Zg wird ein m Fig. 3 dargesteUter 
Fin^l berrachte?hi UR*3 tede? Empol durch eine dicke durchgezogene Unie und the Fornv vonspaler 
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nur FaHe von linearen Einpolen (z. B. Drahten oder Leitungen) untersucbt werden sollen, sondem aach Falle von 
flachig ausgestaheten Einpolen (z. B. Oberflachenbereichen). 

Ji und J2 st ell en die Ausdehniingsfunktionen der Momentenmethode dar. Die Ausdehinmgsfunktionen be- 
schrieben die Form der Stromverteilung in den Einpolen, Die Ausdehnungsfunktionen unterscheiden sicb 
. abhangig von dem Typ der Momentenmethode; wobei jedoch der folgende Beweis fur alle Momentenmethoden 
Gultigkeit besftzt Als Ausdehnongsfunktion kann so wo hi ein sinusfdnniger oder ein sagezahnformiger als auch 
ein impulsf ormiger Stromverlauf verwendet werden. Zudetn trifft der nachfolgende Beweis auch unabhangig 
von einer linearen oder ebenen Form des Einpols zu. Der Beweis wird nachfolgend bezuglich der "stuckweisen" 
smusformigen Momentenmethode dmxhgeftxhrt Es wird jedoch daranf hingewiesen, daB die Ausdehnungsfunk- 
tionen in Fig. 3 der Emfachheh halber dreieckformig durch gestricheke Unien angedeutet sind, wahrend in der 
folgenden Beschreibung von einem sinusfonnigen Stromverlauf ausgegangen wird. 

Bei der "stuckweisen 0 smusformigen Momentenmethode konnen die Ausdehnungsfunktionen der in Fig. 3 
gezeigten Einpole < 1 > bis <4 > wie fblgt ausgedruckt werden. Dabei beschriebt du d*, ds bzw. cU die Lange 
des Einpols <1 >, <2> f <3> bzw. <4>. 

Elektrischer Stromeinpol < t > : Ji = sink(z— zoVsinkdi, 
Elektrischer Stromeinpol < 2 > : Ji = sink(z2 — z)/stnkd2» 
Elektrischer Stromeinpol <3>:J 2 = sink(t— toVstnkcb, 
Elektrischer Stromeinpol <4 > : J2 «= sinkfe — tysinkdV 

Durch Anwendung dieser Ausdehnungsfunktionen wird zunachst die gegensettige Impedanz Zi 3 zwischen den 
Einpolen < 1 > und <3> so wie die gegensemge Tmperianr zwischen den Einpolen < 1 > und <4> ermittelt. 
Diese gegenseitigen Impedanzen Z 13 und Zw sind in Fig. 4B durch die entsprechenden Gleichungen dargestellt 

Unter BerOcksichtigung der Tatsache, daB die Entf ernung zwischen Einpolen durch 

r « (z 2 + t 2 - 2ztcosO + h 2 ) 172 

ausgedruckt werden kann, ergeben sich fur die Koefnzienten dieser gegenseitigen Impedanzen Z13 und Zm die in 
Fig. 5A dargestellten Gleichungen. Es ist zu beachten, daB gilt a = cu/4it (wobei a konstant ist). 

Des weheren kann die Entf ernung r zwischen Einpolen folgendermaSen naherungsweise ausgedruckt werden, 
wobei ro die Entfernung zwischen dem Endpunkt des einen Einpols und dem Endpunkt des anderen Einpols 
angibt: 

r - [r 0 2 + (r 2 - ro 2 )] 172 = (r 0 2 + A) 172 - ro(l 4- A/ro 2 ) 172 

- ro(l + A2r 0 2 — A 2 /8r 0 * + ...) 

= r 0 + A/2r 0 — A 2 /8r 0 3 + . . . = r 0 + d 

(dabeigihA - (r 2 - ro 2 ) und d = A/2r 0 - A 2 /8r 0 3 + . . . und r 0 - (z, 2 + ti 2 - 2z 1 tiCOsO + h 2 ) 172 ). 

Dementsprechend ergeben sich fur die gegenseitigen Impedanzen Z13 und Z M die in Fig. 5B dargestelhen 
Ausdrucke. Dabei gilt Ai = sinkdi x sinkdV 

Urn diese Gleichungen vereinfacht ausdrucken zu konnen, wird u = z — z^v« t - to und w « — t + t2 
gesetzt Des weheren wird die Potynomnaherung sinku = ku — (ku?)/6 und cosku = 1 — Qaif/2 durchgefuhrt 
und angenommen, da6 der Einpol kurz ist so dafi sich fur den ReaheQ Ri der gegenseitigen Impedanz Z13 die in 
Fig. 6 dargesteQte Gleichung und schfieBbch das in Fig. 7 gezeigte Ergebnis ergibt. Es ist zu beachten, daB dabei 
der Anteil e noch nicht berucksichtig ist. 

Werden die Koefnzienten der Terme von k 4 usw. durch P| bis Ps abgekurzt, ergibt sich folgende Gleichung: 

Ri = (0/AiICPik 4 - Pjk 6 + Psk 8 - P4k t0 )cosOi - Psk 2 + Pek* - P 7 k« + Psk*} 

Wird des wetteren naherungsweise "At = sinkdi x sinkdi = did^k 2 " gesetzt und der zunachst unterschlagene 
Anteil e~ 9°° eingesetzt, so ergibt sich fur den Realteil Ri der gegenseitigen Impedanz Zu folgender Ausdruck: 

Ri = [ae-^didak^^ik^-Pjk 6 -f Psk*- P 4 k 10 )cosO l - Psk 2 + Psk 4 - Pyk 6 + Psk 8 ] 
= [ae-^didsHEtPik? - P 2 k* + Pjk 6 - P 4 k 8 )cos<D 1 - P 5 + Psk 2 - P 7 k* + Psk 6 ! 

Der Ixnaginarteil Ij der gegenseitigen Impedanz Z13 wird hingegen durch die in Fig. 8 dargestellte Gleichung 
sowie das in Fig. 9 gezeigte Endergebnis ausgedruckt, wobei wiederum zunachst der Anteil e~^° noch nicht 
berucksichtigt ist. 

Werden die Koefnzienten der Terme von k 3 usw. durch Qi bis Q 8 abgekurzt, ergibt sich folgende Gleichung: 

Ii - (o/AxXWik 3 - 02k 5 + Qjk 7 - Q**)cosOi - Q 5 k + Qek 3 - Qrk 5 + Qgk 7 ]. 

" Wird des weheren naherungsweise *A t = sinkdi x sinkd 2 = d^k 2 * gesetzt und der zunachst unterschlagene 
Anteil e~ j1cr0 eingesetzt, so ergibt sich fur den ImaginarteU It der gegenseitigen Impedanz Zu folgender 
Ausdruck: 

U = [ae-i krt /(d 1 d 2 k 2 )][(Qik 3 - Chk 5 + Q 3 k 7 -Q 4 k9)cosOi - Q 5 k + Qsk 5 - Qrk 5 + Qsk 7 ] 
» [ae-^/CdidaHKQik - Q 2 k 3 4- Qsk 5 - Q+kOcosOt - Qs/k -f Qek - Qrk 3 + Qsk 5 ]. 
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Dementsprechend ergibt sich fur die gegenseitige Impedanz Z«: 

= [M-^/WdhHBPik? - Pzk 4 + P 3 k 6 - P4k«)cosOi - Ps + Psk 2 - Pzk 4 + P«k*] 
+ Jloe-^dids)] [(Qik - Q 2 k 3 + Qsk 5 - Q^oosOH - Qs/k + Qek - Qzk 3 + Qsk 5 } 

Auf ahnliche Weise kann die gegenseitige Impedanz Zu ermittelt werden, wobei sich hierfur die nachfolgende 
Gleichung ergibt, falls die Termevon k 3 etc durch Ri bis Rs ausgedrucktsind: 

2w « [ae-^/CdAHHRi 2 - Rjt 4 + R^k 6 - R^k^cosOh + Rs - R«k 2 + R/^-Rsk 6 ] 
+ j [ae-^dtcU)] [(S 2 k - Ssk 3 + Ssk 5 - S 4 k 7 )cos<I> 1 + Ss/k - Ssk + Sjk 3 - Sek 5 ]. 

Somit kann die gegenseitige Impedanz (Zi 3 + Zm) als Polynom der Welleozahl wie folgt dargestellt werden, 
wobei die Koeffizienten der Terme von k 3 etc. durch C\ bis C9 ersetzt worden sind: 

Zi 3 + Zi4 « e-* 0 [(Co + Qk 2 + C 2 k* + Csk 6 + C* 8 + . . .) + j(Csk- 1 + Csk + C T k* + Cgk 5 + C*k 7 + . . .)} 

Wird mit Z23 die gegenseitige Impedanz zwischen den Einpolen <2> und <3> and mit Z24 die gegenseitige 
Impedanz zwischen den Einpolen <2> und <4> beschrieben, kann auf ahnliche Weise die gegenseitige 
Impedanz (Z23 + Z24) ermittelt werden, so daB insgesamt die gegenseitige Impedanz (Z13 4- Z i4 + Z23 + Z 24 ) als 
Polynom der WeUenzahl k durch eine zu der obigen Formel ahnliche Gleichung dargestellt werden kann. 

Die gegenseitige Impedanz Zjj kann demnach allgemein wie folgt dargestellt werden: 

Zy - e-i^gao + atk 2 +a 2 k* + ask 6 + ink 8 + ...) + Kbok- 1 4- b,k + b* 3 + bak 5 + l*k 7 4 ...)J 

Es kann bewiesen werden, daB auch die gegenseitige Impfdanz Yg analog zu der gegensem'gen Admittanz 
allgemein durch ein Polynom der Frequenz dargestellt werden kann. Die Naherungsgleichung fur die gegenseiti- 
ge Admittanz Yij entspricht genau derjenigen fur die gegenseitige Impedanz Zy Aus diesem Grande wird auf den 
Nachweis der entsprechenden Naherungsgleichung verzichtet 

Die gegenseitige Impedanz Yij kann demnach ebenso allgemein wie folgt dargestellt werden: 

Yy « c-i^Rao + aik 2 + a 2 k« + 33k* + . . .) + jCbok' 1 + b,k + t^k 3 + bjk 5 -f 04k 7 4- ."..)]• 

Nachfolgend wird die DarsteUung der Wechserwirknng Bij durch eine Polynomfunktion der Frequenz erlau- 
tert 

Zu diesem Zweck wird der in Fig. 10 dargesteUte Einpol betrachtet In dieser Figur ist der Einpol durch erne 
dicke durchgezogene Linie dargestellt, wihrend die strichpunktierte Linie analog zu Fig. 3 die Form der 
Ausdehnungsfunktion angibt 

Es wird das aufgrund der Stromquelle hervorgerufene Magnetfeld untersucht, wobei die Emheitsvektoren wie 
in Fig. 1 1 A gezeigt definiert ist 

Gernifi der stuckweisen sinusformigen Momentenmethode konnen die Ausdehnungsftinktionen der in Fig. 10 
gezeigten Einpoie < 1 > bis <4> folgendermaBen beschrieben werden, wobei di, da dj bzw. d* die Lange des 
Einpols < 1 >, <2>, <3> bzw. <4> bezeichnet: 

Elektrischer Stromeinpol < 1 > : Ji = smk(z— zo)/sinkdi. 
Hektrischer Stromeinpol < 2 > : J 2 = sinkfe -z)/sinkd2. 
Magnetischer Stromeinpol <3>:M3 = sink(t— to)/sinkd3. 
Magnetischer Stromeinpol <4>:M4 = sink(t 2 — t)/sinkcU- 

Dabei ist die Stromquelle ausschlieBlich auf der z-Achse angeordnet, das elektromagnetische Feld ist (unab- 
hangig von der O-Koordinate) kreissymmetrisch ausgestaltet und es tritt nur ein H^-Antefl auf. Dieser HotAd- 
tefl ist durch die in Fig. 1 IB gezeigte Gleichung definiert 

Bei der in Fig. 10 gezeigten Anordnung betragt die tangenuale Komponente des Magnetfeldes des Einpols 
<3> Ht = — (h/p) x H(*sind>, so daB die Wechselwirkung Bij durch die in Eg. 11C dargesteUte Gleichung 
beschrieben werden kann. m . 

Werden die in Fig. 12 dargesteuten Bedingungen in diese Gleichung eingesetzt, ergibt sich die in Fig. 13A 
dargesteUte Gleichung fur die Wechselwirkung Bij. 

Des weiteren kann , wie bereits zuvor beschrieben worden ist, die Entfernung r zwischen den Einpolen 
ausgedruckt werden durch: 

r = (z 2 + 1 2 - 2ztcosO + h 2 ) 1 ^ - r 0 + d. 

Fur die Wechselwirkung Bij ergibt sich demnach die in Fig. 13B gezeigte Gleichung. 

Ausgehend von dieser Gleichung kann die Wechselwirkung Bn zwischen den Einpolen < 1 > und <3 > so wie 
die Wechselwirkung B M zwischen den Einpolen < 1 > und <4> gemaBFig. 14 beschrieben werden. 

Urn diese Gleichungen vereinfacht ausdrucken zu konnen, wird u-z-zo,v = t- to und w = — t + t 2 
gesetzt Des weiteren wird die Porynoranaherung sinku = ku — (ku 3 )^ und cosku « 1 - (kuf/2 durcfagefuhrt 
und angenommen, daB der Einpol kurz ist, so daB sich fur den zweiten Term auf der rechten Sehe der 
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Wechsdwirkurigsgleichung der in Fig. 15 gezeigte Ausdruck und schlieSIich das in Fig. 16 gezeigte Ergebnis 
ergibt 

Werden die Koefftzienten der Terme von k 2 usw. durch Pi bis Pa abgekurzt, kann die in Fig. 16 gezeigte 
Gieichung folgendennaBen umgescbrieben werden: 

Pik 2 + P 2 k* + P3k« + P^k 8 + jCPsk 3 + PeK 5 + rVk 7 + Pgk 9 )- 

Auf ahnEche Weise kann aach der erste Term auf der rechten Seite der die Wechseiwirkung B13 beschretben- 
den Gieichung in die in Fig* 17 gezeigte Gieichung und schliefllich in den in Fig. 18 gezeigten Ausdruck 
umgewandeft werden. 

. Werden die Koeffizienten der Terme von k 2 usw. durch Qi bis Qs abgekurzt, kann die in Fig. 18 gezeigte 
Gieichung foigendermafien umgeschrieben werden: 

Chk 4 + 02k 6 + Oak 8 + Q4k»> + flQsk 3 + Qek 5 + Q* 7 + Qsk 9 ). 

Mit diesen Abkurzungen kann die Wechsebrirkung B13 schlieBfich f olgendermaSen ausgedruckt werden: 

B„ « I>isinO,/(4jismkd| x sinkdsfr-* 10 x [Ptk 2 + P 2 k* + Pak 6 + P4IC* + j(Psfc* + Pek 5 + P 7 k 7 + Psk 9 ) 
+ Qik 4 + Qz*? + Oik 8 + Q* 10 + KQ5IC 3 + Q«k* + Q* 7 + Qsk 9 )} 

Wird des weheren naherungsweise sinkdi x sinkd3 = didjk 2 eingesetzt, so ergibt sich fur die Wechseiwirkung 
B13 folgende Gieichung: 

B l3 = PusinOj^mlids)^-* 10 x {P, + CP 2 + Qi)k 2 + (P 3 + Q^k 4 + (P 4 + Qajk* + Q^k 8 

+ I(Ps + Q 5 )k + (P 6 + Qejk 3 + (P 7 + Qzjk 5 + (P* + Qs* 7 ! 

= e-^fRt + Rak 2 + 1^ + 1*4^ + Rsk 8 + j(R«k + Rjk 3 + R«k 5 + Rsk 7 )} 

Auf ahnliche Weise kann die Wechseiwirkung Bu ermhtelt werden, wobei sich hierfur die nachfolgende 
Gieichung ergibt: 

Bm - e"** [Si + Sjk 2 + S^k 4 + S4k 6 + Ssk* + j(S«k + 87^ + 58^ + S&)1 

Somit kann die Wechseiwirkung (B13 + Bt4) als Pc4ynom der Wellenzahl wie foigt dargestelh werden: 
B« + B14 - e-^Cd + C 2 k 2 + Cjk 4 + 04k* + Csk« + j(Cek + Cjk 3 + Cgk 5 + C^k 7 )]. 

Wtrd mit B23 die Wechseiwirkung zwischen den Einpolen <2> und <3> und mit B** die Wechseiwirkung 
zwischen den Einpolen <2> und <4> beschrieben, kann auf ahnliche Weise die Wechseiwirkung (B23 + Bm) 
ermhtelt werden, so daS insgesamt die Wechseiwirkung (B 13 + Bu + B23 + B24) als Pofynom der Wellenzahl k 
durch eine zu der obigen Formel ahnliche Gieichung dargesteDt werdenkann. 

Die Wechseiwirkung Bij kann demnach aOgemein wie folgt dargestelh werden: 

Bij - e-^Kao + atf 2 + a 2 f* + asf* + a+f 8 + . . .) + Kbof 1 + bif 3 + fcf 5 + bjf 7 + ttf 9 + ...)} 

Nacfafolgend wird die vorliegende Erfmrinng detaiDiert anhand der von der erfindungsgemaBen Simulauons- 
vorrichtung 1 chuxhgefuhrten Verfahrensschritte erlautert Es wird darauf hingewiesen, daB der Einfachheit 
halber die Simulation bei dem nachfolgend beschriebenen Ausfuhrungsbeispiel ohne Berucksichtigung der 
gegenseitigen Admittanz und der Wechseiwirkung durchgefuhrt wird. 

Fig. 19 zeigt ein Ausfuhrungsbeispiel der von der erfmdungsgemaBen Simulationsvorrichtung 1 durchgefuhr- 
ten Verfahrensschritte. Dabei ist mit 100 eine Verwaltungsdatei bezeichnet, welch e die Konfigurationsinforma- 
tionen der zu simulierenden und in Elemente unterteUten elektronischen Vorrichtung verwahet 

Wie aus diesem FhiBdiagramm erstchtUch ist, wird in Schritt 1 (ST1) zunachst eine "Eingabe einer Wellenquel- 
le", in Schritt 2 (ST2) eine "Berechnung einer Naherungsgieichung*, in Schritt 3 (ST3) eine "Berechnung eines 
Oberwellenstromes", in Schritt 4 (ST4) eine "Berechnung eines Stromes*, in Schritt 5 (ST5) eine "Berechnung 
einer Spannung", in Schritt 6 (ST6) eine "Berechnung eines dektromagnetischen Feldes", in Schritt 7 (ST7) eine 
berechnung einer zwischen zwei Leitern auftretenderj Spannung* und in Schritt 8 (ST8) eine "Berechnung der 
Kopplung einer Obertragungslettung" diirchgefuhrL 

Fig. 20 zeigt ein Ausfuhrungsbeispiel der "Eingabe einer WellenqueUe" gemaB des in Fif. 19 gezeigten 
Schrittes ST1. 

Wie in diesem FluBdiagramm dargestelh ist, werden wahrend des Schrittes "Eingabe einer Wellenquefle" hn 
Zehbereich vorliegende Wellenquellendaten aus der Verwaltungsdatei 100 ausgelesen. 

Wie in Fig. 21 gezeigt ist, werden ausgelesen Daten unterteilt und fouriertransformiert, urn auf diese Weise 
das Frequenzspektrum der WellenqueUe zu erhahen. Wird zu diesem Zeitpunkt der Gleichanteil dieses Fre- 
quenzspektrums bendtigt, kann nach dem Ohm'schen Gesetz dieser Gleichanteil erhahen werden. 

Zum Teste n einer elektronischen Vorrichtung kann ein Hochspannungsimpuls oder ein hoher Stromimpuls 
angelegt und dessen Auswirkungen untersucht werden. Zur Simulation dieses Vorganges muB als WellenqueUe 
ein mit einer Anfangsspannung geladener Kondensator oder eine mit einem Anfangsstrom versorgte Spule 
angenommen werden. 
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Fur den Fall, daB die WeUenqueUe durch einen mit einer Anfangsspannung geladenen itondensator geb^t 
wirTl^ dlSeWcUenquelle in die in Kg. Z2A gezeigte Ersatzschaltung urngewandelt werden. Da ffir die 
^erXg der FourierLnsformatioii Periodizitat erforderlich ist.wird.wfe .n Fig. gezejgt ist, erne von 
^ E^a&n^lSe stufenfonnig verlaufende Spannung der Etafachhert ^^pannungam- 
Dutemit langen Perioden 1 1 und t2 behandelt Durch die Anfangsspannung des Kondensators wird imallgemei- 
SeTeto gXpfte^hwingung hervorgerufeu. wobei die Zeh U. zum Au^w^en d,es« g^pft^ 
SSonvorhemgbar isTso daB den Zeitspannen tl und 12 em Wert zugewiesen wrrd. der groBer als diese 

die* Wdleuqueue in die in Fig. 23 A gezeigte Ersatzschaltung urngewandelt werdea Da w.ederum to die 
Anwendung der Fouriertransfomation Periodizitat erforderlich ist, wird, wie m K*23B gezagt ist, em yon 
erzeugter stufenfonnig vcrlauf ender Strom der Enfachhert halber ak Strommipulse rnit 
^ip^riodeTtluSl t2 behandek. Durch den Anf angstrom der Spule wnd im aUgememen erne gedampfte 
SServorgerufen, wobei die Zeit bis zum Ausschwingen dieser gedampften OsziBauon vorhersagbar 
^^Hen Zeitspannen ti und t2 ein Wert zugewiesen wird. der groBer afe diese zum Ausschwmgen der 

^A^^WeS wS'w^nd des Schrittes -Eingabe einer Weflenquelle- die im Zehbereich voruegenden 
Wef r^eUe^n Sfder Verwaltungsdatei 100 ausgelesen, unterteflt und fouriertnmsfornuert, um somit das 
F^aucn>»£k^dc? entsprechendefwellenquene zu erhalten. AnschlieBend wird. wje beretts erwahnt 
wmteufetmitA^ Schritt ^rechmmg einer Naherungsgleichung" (ST2 in Fig. 19)fortgefahren._^ 
^fze? eta A^rungsbeispW fur dec Ablanf des in F.g. 19 gezeigten Schnttes -Berechnung emer 

N teT1e^Bdiagnunm dargesteUt ist. wird wihrend der -Berechnung einer 
znnach^in Schritt 1 (STl)ein Frequenzbereich aus den P<»^ s 1 f^^ en *^ h ££*2££5 
S einer WeUenqueUe- erhaltenen Frequenzspektrums der WeUenqueUe ausgewahlt bzw . speznmert 1st 
tone GeSceh erforderlich, wW der gesamte Frequenzbereich ausgewahlt SoU hmgegeu pnmar 
dTsh^atior^eAum werden, kSunen die OberweUenanteUe der Frequenz entfemt ^^* e 
Sn«^Z*len Frequenzbandes Hegenden Frequenzen ausgewahlt bzw. fur die Simulation verwendet werden. 
S TathSend 3S^tt^^AbS*eq^enaus dem "Sc^^^f^"^^ 
ausK^ablt So konnen beispielsweise fflnf Abtastfrequenzen f,„ U fa U bzw. fa ausgewahlt werdetL Die 
XSSue^werden^ dem in dem Schritt MST!) Testgelegten 

Hnreh siedieser Frequenzbereich grob unterteik wird. Fur erne hohe Genaingkeit muBdie Anzanlder Abtes»e- 
H^tSSt we^I^en* die Anzahl der Abtastfrequenzen zu reduzieren, falls pnmar erne Verkur- 

^E£SK2C± - den bisher no* ni^verarbe^ Ab^enzen m Scbritta 
(Sn)wWruTdie ausgewahlte Abtastfrequenz in Schritt 4 (ST4) die geger*e.^tap«lanz Zij & J T*^- J*? 
HOfeetaS gen.uen B^echnungsverfahrens berechnet, und der berechnete We* wu£m emer Arbertsdate. 200 
^eicher£ AnschlieBend wird wieder zu dem Anfang des Schrittes 3 (ST3) zuruckgekehrt 
^^chdem jedetetgdegte Abtastfrequenz verarbehet worden ist. wird in Schntt 3 (ST3) erkannV daB die 
aer gegeieitigen Impedanz Zg (i. j = 1 ... M) fur alle Abtastfrequenzen ar^esc^en i* und « 
^S^ritt 5(ST5) foriefahren, in dem die folgende NaherungsgleichunB fur <ke gegensemge Impedanz Z., 
bestimmt wird: . , 

fwnhei f die Frecmenz, k die Wellenzahl bzw. r 0 den Abstand zwischen zwei Hementen i und j beschreibt> Die m 
tteS -Ki die entsprechende Werte der ^^^^ ^ 

warden in diese Gleichung eingesetzt, um auf diese Weise Simu^auoiisgleichungen 

Tund bo bis b L als Unbe^nnte zu erhalten. Diese Simulatioiisgieichungen werden gelos^um dieKoeffiz^eaten 
LbTaLund^bisbLZuenmtteln.DieseK^ 

Wurde^beUpidsweise die gegenseitigeii Impedanzen 25j sl bis Zij* zwischen den Elementen i und j fur funf 
Abtastfrequenzen fsl bis fs5 folgendermaBen ermitteh: 

Zijsi — cd^i + jpijsi, 
. Zij$2 = ai]j2 + jpijja 

Zijs3 = Ctijs3 + jP»j$3, 
ZjjM = cuj$4 + j#js4, 
Ziis5 = aij s s+jPiM 

so ergeben sich die folgenden Simulationsgleichungen: 

ai U - exp[(-j(Wc)ro][ao + aiW + + arf*aj + 

ou U = exp[(-j(Wc)ro][ao + a,tf + + ^ + ^ 

oiL - exp[(- j(2rf s4/c)ro] [ao + a,f* 2 + a*W + a^« + arf^ 

und 
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Pijsi - cxp[(-j(2nf,i/c)ro][bo + M $ i 2 + bzhi* + Wsi 6 + b4n* , 
Pij*2 - exp[(-j(2Kf 52 /c)ro][bo + biftj 2 + bate 4 4- Wsz 6 + b4if , 
Pij*3 = exrf(-j(2«fi3/c)ro][bo + bif,^ + bzfss* + bate 6 + b 4 f s3 !r , 
pijs* * exrfC-jtMs^cJroJEbo + bifs* 2 + bzk* 4 + Ws4 6 + b4s* r , 
PijsS = eip[(-j(2jxf S 5/c)ro][bo + W*? + + WjS 5 + b^sS* . 

Durch Losen dieser Simnlatio nsgleichungen konnen die Koeffizienten ao bis as and bo bis b* ennittelt werden, 
so daB sich die folgende Nahenmgsgleichimg fur die zwischen einem Element i und einem Element j auftretende 
gegenseitigeimpedanz Zij ergibt: 

Zg = e-^Rao + aiP + a 2 f* + asf 6 + ) + Xbof 1 + oif + bjf 3 + hjf 5 + b^f 7 )]. 

Es ist zubeachten, dafi diese Nahenmgsgleichung nnahhangig fur jede Frequenz ausgewertet werden kann, so 
daS eine wehere Erhdhung der Verarbeitungsgeschwindigkeit durcb die Verwendung eines parallelen Compu- 
ters erzielbar ist 

Nach der Ermhtlung der geeigneten Nahenmgsgleichung fur die zwischen Elementen i and j auftretende 
gegenseitige TmpoHati-r Zij (i, j «* 1 ... M) durch den Schritt "Berechnung einer Naherungsgleichung" wird mit 
dem Schritt "Berechnung eines Oberweflenstromes" (ST3 in Fig. 19) fortgef ahren. 

Fig* 25 zeigt ein Ausfuhnmgsbeispiel eines FhiSdiagrammes fur den Schritt "Berechnung eines Oberwellen- 
stromes". 

Wie aus diesem FhiBdiagranun erstchtlich ist, wird in Schritt 1 (ST1) der "Berechnung eines Oberwellenstro- 
mes" zunachst eine noch nicht verarbeitete Frequenz aus den positiven Seitenfrequenzen des in dem Schritt 
"Eingabe einer Wellenquelle" erhahenen positiven Seitenfrequenzspektruins ausgewahlL 

Nachf olgend wird in Schritt 2 (ST2) die ausgewahlte Frequenz in die Naherungsgleichung fur die gegenseitige 
Tmpffdanz Z«j emgesetzt, welche in dem Schritt "Berechnung einer Naherungsgleichung" ermitteh worden ist, urn 
die zwischen einem Element i und einem Element j auftretende gegenseitige Tmpedan7 Zij (i, j = 1 ... M) zu 
berechnen. Die Berechnung dieser gegensehigen Impedanz ist ein emfacher Einsetzvorgang, so daB die Berech- 
nung mit einer auBerst hohen Geschwmdigkett durchgefuhrt werden kann. 

AnschfieBend werden in Schritt 3 (ST3) die in Fig. 26 gezeigten Simulationsgleichungen gemaB der Momen- 
tenmethode geldst, indem die in Schritt 2 (ST2) berechnete gegenseitige Impedanz Zij (i, j « 1 . . . M) and das 
Frequenzspektrum V(e>) (co = 2nfq) der WeDenqueDe fur die in Schritt 1 (ST1) ausgewahlte Frequenz fq 
ausgewertet werden, urn die in jedem Element bei der in Schritt 1 (ST1) ausgewahhen Frequenz f q flieBenden 
Strdme Ii(fq) bis lM(fq) zu ermitteln. Dieser Vorgang wird wiederholt, bis samtfiche Frequenzen in Schritt 1 (ST1) 
ausgewahlt worden sind. Auf diese Weise werden fur samtliche Frequenzen f q (q = 1 . . . N) die in jedem Element 
flieBenden Strdme Ii(fq) bis InM(£q) ermitteh. 

Dabei wird bei den in Fig. 26 gezeigten Simulationsgteichungen der M omentenmethode angenommen, daB 
die Weflenquelle V(<a) in dem Element m vorhanden ist. Es sei darauf hingewiesen, daB eine we it ere Verarbei- 
tungsbeschlennigung erzieh werden kann, wenn bei der Losung der Simulationsgleicbungen der Momentenme- 
thode beispielsweise die (in der Japanischen Patentanme Idung 7-342695 offenbarte) LU-Aufgliederungsmethode 
verwendet wird. 

Das in dem Schritt "Eingabe einer Weflenquelle" erhaltene Frequenzspektrum weist, wie aus Fig. 21 erstcht- 
lich ist, auch ein negatives Seitenfrequenzband auf. Wie oben beschrieben worden ist, werden durch die Schritte 
"Berechnung einer Naherungsgleichung* und "Berechnung eines Oberwellenstromes" lediglich die von einer 
Wellenquelle mit positiven Seitenfrequenzen indozierten Strdme berechnet. Daher mussen auch die von einer 
Weflenquelle mh negativen Settenfrequenzen induzierten Strdme ermhtelt werden. 

Diese Berechnung wird durchgefuhrt, indem die Tatsache ausgenutzt wird, daB zwischen den in jedem 
Element m flieBenden Stromen Im(co) und Im(— <a) ein konjugiert-komplexer Zusammenhang besteht, da zwi- 
schen Zjj(6)) und Zij( — <a) der folgende konjugiert-komplexe Zusammenhang vorhanden ist: 

Zij(e>) « e-^Rao + aif 2 + a 2 f* + ajf 6 + a^f*) + jCbof 1 + btf + b^P + baf 5 + b^f 7 )] 

Zij(-<a) = [(ao + aif 2 + a 2 f* + aaf 6 + a^f 8 ) - Xbof 1 + bif + bzf 3 + bsf 5 + W*)]. 

Zwischen den Frequenzspektren V(a>) und V( — <a) der Wellenquelle besteht ebenfafls ein konjugiert-komple- 
xer Zusammenhang. Der durch eine Weflenquelle mit negativen Seitenfrequenzen hervorgerufene Strom 
I m (_e)) kann demnach durch Umkehren des Vorzetchens des Imaginirteils des von der Wellenquelle mit den 
positiven Seitenfrequenzen hervorgerufenen Stromes I m (co) ennittelt werden. 

Es ist zu beachten, dafi diese "Berechnung eines Oberwellenstromes" unabhangig fur jede Frequenz mogfich 
ist, so dafi dieses Verf ahren durch die Verwendung eines parailelen Computers weiter bescfaleunigt werden kann. 

Nach Berechnen des in jedem Element bei entsprechendeo Frequenzen flieBenden Stromes mit Hilfe des 
Schritt es "Berechnung eines Oberwellenstromes" wird das Verf ahren mh dem Schritt "Berechnung eines Stro- 
mes* (ST4 in Fig. 1 9) f ortgesetzt 

Fig. 27 zeigt ein Ausfuhrungsbeispiel eines FluBdiagramms fur diesen Schritt "Berechnung eines Stromes". 

Wie in diesem FluBdiagramm dargestellt, werden zunachst im Laufe des Schrittes "Berechnung eines Stro- 
mes", nachdem durch einen Benutzer derjenige Ptinki (bzw. dasjenige Element) bestimmt worden ist, fur den der 
Strom berechnet werden soli (ST1), die bei entsprechenden Frequenzen durch diesen bestimmten Punkt flieBen- 
den Strome "Ip(fi) bis Ip(fN), lp(0) und Ip(— fi) bis Ip(— tuff aufgegriffen bzw. ausgelesen und anschlteBend die 
inverse Fourier-Transformation durchgefuhrt (ST2), um den in dem bestimmten Punkt p flieBenden Strom im 
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Zeitbereich zu berechnen. -. . . . . w . t>,,„u + 

Nachdem mit Hilfe dieses Verfahrcnsschrittes -Berechnung ernes Stromes" em in emem bestimmten Punkt 

flieBender Strom im Zeitbereich berechnet worden ist, wird mit dem Schritt berechnung einer Spannung (vgi 

STSinFig. 19)fortgefahreiL . . . 

Re. 28 zeigt ein Beispiel fur die Durchfuhrung der berechnung emer Spannung'. 

Wie in diesem FhiBdiagramm dargestellt ist, werden im Laufe des Schrittcs berechnung emer Spannung* , 
nachdem durch einen Benutzer derjenige Punkt (bzw. dasjenige Element) bestimmt worden *t,far den die 
Spannung berechnet werden soil (ST1X zunachst die bei entsprechenden Frequenzerj idurch diesen bertunmten 
pSflfeBenden Strome Tpfc) ^ l^l «0) und Ip(-f,) bis IpC-f^ aufgegrrfto tew. ausgeIesen£T2) mid 
SeBendmeseWertemitd^ 

Punktes p for die einzelnen Frequenzen multipliztert, um die in dem besummten Punkt p bei den einzelnen 
FrequenienauftretendeSpann^^ 

Nachfolgend werden diese somit fur die einzelnen Frequcnzen ermittelten Spannungen 'V^U) i bis 
Vp (o) und Vp(-fi) bis V^-faF einer inversen Fourier-Transformation (ST4) unterzogen, um die in dem 
besummten Punkt p auftretende Spannung im Zeitbereich zu berechnen. . ^ . 

Nachdem mh Hflfe des Schrittes -Berechnung einer Spannung* die in emem bestmunten Punkt auftretende 
Spannung im Zeitbereich ermittelt worden ist, wird das Verf ahren mit dem Schritt -Berechnung ernes elektroma- 
gnetischenFeldes'fortgesetztCvgLSTSinFifrl^ _ , 

' Fig. 29 zeigt ein Ausfuhnmgsbeispiel fur ein RuBdiagramm fur den Verfahrensschntt Berechnung ernes 

^r^eTwWz^Erfeute^ der Fall betrachtet, daB die Intensitat eines elektrischen Feldes berechnet werden 
solL Wie in diesem FraBdiagramm dargestellt ist, wird zunachst in dem Verfahreiisschrm^rechnung ernes 
elektromagnetischen FeWes" durch einen Benutzer ein Untersuchungspunkt p festgelegt (ST1). AnschheBend 
wird aufgrund bekannter theoretischer Formem fur elektromagnetische Felder das an dem Untersiicnungspunkt 
ry auftretende und durch die in jedem Element bei den einzelnen Frequcnzen flieBenden Strome hervorgerufene 
elektrische Feld ermittelt, umauf diese Weise die in dem Untersuchungspunkt p bei den besummten Frequenzen 
erzeugtenelektrische Felder bp(fi)bis Etf*l E^O)imd^-f } )bisEp( r W^^k , 

Nachfolgend wird auf die fur die einzelnen Frequenzen ermittelten elektnschen Felder T^U) bis E^hi} E£0) 
und EoT-fi) bis EJ-tuT die inverse Fourier-Transformation angewendet, um auf diese Weise das in dem 
^r^ch^unkt p erzeugte elektrische Feld im Zeitbereich zu berechnen (ST4). Auf ahnhch Weise wird 

^^B^^d^d^^^etisd^ Feldes mh HOfe des VerfahrensschrittesJ^ ernes eiek- 

tromagnetischen Feldes-, wird mh dem Verfahrensschntt berechnung einer zwischen Lertem anhegenden 
Spannung* (ST7 in Kg. 19) fortgefahren. , _ _ . , ^ 

Bei diSem Verfahrensscfaritt berechnung einer zwischen Leitern anliegenden Spannung- kann nach Festle- 
gung einer zu verarbertenden SteUe p zwischen Leitern durch einen Benutzer die zwischen den Lertem bei emer 
Snmten Kreisfrequenz <* auftretende Spannung VrfeO durch .folgende Fomd J^f^ w ^ h ^ 
der bei einer besummten Kreisfrequenz ffl in einem Element n flieBende Strom nut Up) mid die zwischen der 
Stelle p der Leiter und dem Element n auftretende gegenseitige Impedanz nut Zp^o) bezeicnnet ist: 

Fonnel 1: V p (m) = -f) «<d) ZJp). e 
«=i 

Auf diese Weise kann die bei besummten Frequenzen zwischen den Leitern anliegende Spannung "Vp(f i) bis 

V ^c^Be^ d wu4"au] d^ri dLn ^e^mte^Fr^quenzen zwischen den Leitern anliegenden Spannungen 
*Vt/fi)bis VpffK), V„(0)und Vp(-fi) bis V^-fxY die inverse Fourier-Transformauon angewendet, um somit die 
zwischen den Leitern an der Stelle p auftretende Spannung im Zertbereicfa zu erhalten. 

Zur Eriauterung der obigen Formel 1 wird Rg. 30 betrachtet Ist zwischendie mrt pi und p2 benfao^ 
Leiter ein Widersiand R eingefugt, ergibt sich die in Fig. 31 A dargestellt* Gleichung 

gung, daB das elektrische Feld um jedeh Leiter Null ist Daraus kann die m Fig. 31B gezeigte f^chung fur den 
zwischen den Leitern flieBenden Strom Ip abgeleitet werden, so daB sici wiedennn fur die zwischen den Lertern 
arikgende Spannung Vp die in Fig. 31C gezeigte Formel ergibt Tatsachhch fiieBt jedoch zwischen den Lertem 
Strom, so daB in die in Fig. 31C dargesteUte Formel die Werte "R— und Ipi, Ip2 -~o- emgesetzt werden 
konnen. Dadurch ergibt sich die oben angegebene Formel 1. 

Die Formel 1 wird somit aus dem von einem imaginar zwischen die Letter exngefugten Widerstand erzeugten 
Spannungsspektrum abgeleitet, wobei der Widerstand ucendiich groB gewahlt wird 

Nachdem auf diese Weise die zwischen Leitern anliegende Spannung mrt HOfe des Verfahrensschrrttes 
-Berechnung einer zwischen Leitern anliegende Spannung* ermittelt worden ist wird mit dem Verfahrensschntt 
berechnung der Kopplimg emer Obertra^ , _ 

Der Schritt berechnung der Kopplung einer Obertragungsleitung- dient zur Simulation des einer ^Pfangs- 
schaltung zugefuhrten Stromes bzw. einer entsprechenden Spannung, die durch em externes elektnsches Feld in 
einer fiber einer weiten Bodenebene angeordneten Obertragimgsleitung hervorge^enen werde^ 

Dabei sei die Obertragungsleitting wie in Fig. 32 dargestellt uber der weiten Bodenebene an geordne^L Wird 
mit Za die Impedanz ant sendeseitigen Ende. mit L die Uinge der Obertragungsleitung.^t h fi^^^^ 
Obertfagungsleitung, mit p das Obertragungsmafl bzw. der Fortpflanzungskoeffizient der Obertragimgsleitung, 
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mil Zo der WeUenwiderstand der Obertragungsleitung, mit Ex^x, z) der x-Antefl in der Stetle (x, y) der Obertra- 
gungsleitung auf grand des in einem anderen Leiter flieBenden Stromes wirkcnden elektrischen Feldes und mit 
Ez*(x, z) der z-Anteil des in der Stelle (x, y) der Obertraguiigsleitung aufgnmd des in einem anderen Leher 
flieBenden Stromes wirkenden elektrischen Feldes bezeichnet, ergibt sich fur K(x) folgender Ausdruck: 

5 Fonnel2:K(x) = Er i (x t h)-ExXx f 0). 

FurDgilt: 

10 Fonnel3:D=(2 0 Zd+ZoZr)cospL+j(Zo 2 + ZdZr)sinPL fc 

Ez* and K(x) konnen aufgnmd der auf anderen Leitem auftretenden StromverteUung sowie der mit Hilfe der 
Momentenmethode berechneten gegensehigen Impedanz durch die in Fig. 34 gezeigte Fonnel berechnet wer- 
den, wenn die Obertragungsleitung wie in Fig. 33 gezeigt durch die smusfdnnige Ausdehnimgsfunktion Ji usw. 
15 ahnlich zu der Momentenmethode unterteih bzw. beschrieben wird. Dabei bezeichnet Zom die zwischen der 
oberen Leitung der Obertragungsleitung und einem anderen Leiter anftretende gegenseitige Impedanz und 
Zom' die zwischen der unteren Leitung der Obertragungsleitung und dem anderen Leiter aoftretende gegenseiti- 
ge Impedanz. 

Es ist auch bekannt, daB bei der Existenz eines derartigen externen elektrischen Feldes Ex 1 und Ez 1 der an dem . 
2o sendeseitigen Ende der Obertragungsleitung auftretende Strom Id(e>) durch die in Fig. 35 dargestellte Formel 
von "C. D. Taylor" ermhtelt werden kann. Der Term A dieser Gleichung ist in Fig. 36A aufgelost dargestelh, so 
daB sich fur den am sendeseitigen Ende der Obertragungsleitung auftretenden Strom Id(a>) der in Fig. 36B 
gezeigte Ausdruck ergibt. . 

Des weiteren ist bekannt, daB bei der Existenz eines externen elektrisches Feldes der am empfangssehigen 
25 Ende der Obertragungsleitung auftretende Strom Ir(o) durch die in Fig. 37 gezeigte Gleichung von "C D. 
Taylor" beschrieben werden kann, wobei sich fur diesen Strom schliefilich der in Fig. 38 dargestellte Ausdruck 
ergibt. 

Es wird darauf hingewiesen, daB diese Formem von "C D. Taylor" in "The response of a terminated two-wire 
transmission line excited by a nonuniform electromagnetic field*, C p. Taylor, R. S. Satterwite und C W. 
30 Harrison, Jr, IEEE Trans. Antennas Propagation, Ap-13, No. 6, Sehen 987—989, Nov. 1967 offenbart sind. 

Somit kann der am sendeseitigen Ende der Obertragungsleitung auftretende Strom Id(6») und der am emp- 
fangsseitigen Ende der Obertragungsleitung auftretende Strom Ir(co) anhand der von "C D. Taylor 9 ermittehen 
Formem bestimmt werden, wenn das externe elektrische Feld Ex 1 und Ez* (bzw. Ez 1 und K(x)) hinsichtlich der 
Obertragungsleitung bekannt ist, und die ermittelten Stromwerte werden mit der Impedanz multipiiziert, um auf 
35 diese Weise die am sendeseitigen Ende der Obertragungsleitung anfiegende Spannung Vd(eo) und die am 
empfangssehigen Ende der Obertragungsleitung anfiegende Spannung Vr(e») zu ernuttem. 

Aufgnmd der von *C. D. Taylor" entwickelten Gletchungen kdnnen daher die am sendeseitigen bzw. emp- 
fangssehigen Ende der Obertragungsleitung bei bestimmtea Frequenzen auftretenden Strome *Ip(fi) bis Ip(fN), 
lp(0) und Ip(— fi) bis l^—hif durch den Verfahrensschritt "Berechnung der Kopplung einer Obertragungsiei- 
40 tung" ermittelt und anschlieBend der inversen Fourier-Transformation unterzogen werden, um die am sendesei- 
tigen Ende bzw. empfangsseitigen Ende der Obertragungslertung auftretenden Strome im Zehbereich zu be- 
rechnen. Des weiteren werden die am sendeseitigen bzw. empfangsseitigen Ende der Obertragungsleitung bei 
. . bestimmten Frequenzen anUegenden Spannungen *Vp(fi) bis Vp(fN), *Vp( 0 ) und Vp(— fi) bis Vp(— fsT Hilfe 
der Gleichungen von "C, D. Taylor* ermittelt und der inversen Fourier-Transformation unterzogen, um die an 
45 dem sendeseitigen bzw. empfangsseitigen Ende der Obertragungsleitung an&egenden Spannongen im Zeitbe- 
reich zu emntteln. 

Durch den Verfahrensschritt "Berechnung der Kopplung einer Obertragungsleitung" kann daher der einer 

Empfangsschaknng zugefuhrte Strom bzw. eine entsprechenden Spannung, die durch ein ext ernes elektrisches 

Feld in einer Obertragungsleitung hervorgeruf enen werden, shuaKert werdea 
so Fig. 39 A und 39B zeigen das Ergebnis der Simulation infolge des Schrittes "Berechnung der Kopplung einer 

Obertragungslehung*. Diese Simulation wurde unter den in Fig. 40 gezeigten Bedingungen dHnxhgefuhrL 
Fig. 39A zeigt den Seodestrom einer angenommenen Sendeantenne, wahrend Fig. 39B den infolge des Sende- 

stroms fiber die Obertragungsleitung flieBenden Strom darstellt. Bei dieser Simulation wurde simutiert, daB der 

in einer von der Sendeantenne lm entfernt angeordneten Obertragungsleitung flieBende Strom mit demselben 
55 Wellentyp bzw. (Mode) und mit einer der Obertragungszeit der Radiowellen entsprechenden Zekverzogenmg 

flieBt, so daB die Gultigkeit des Verfahrensschrittes "Berechnung der Kopplung einer Obertragungsleitung" 

uberpruft werden konnte. 

Wird eine sich zeitlich andernde Wellenquelle vorgegeben, transformiert die erfmdiingsgemaBe Simulations- 
vorrichtung 1 diese Wellenquelle in den Frequenzbereich, ermittelt auBerst schnell die gegenseitige Impedanz 
6o durch Verwendung der Naherungsgleichung, simuliert im Frequenzbereich die Intensitat des elektromagneti- 
schen Feldes usw. mit Hilfe der Momentenmethode und fuhrt hinsichtlich des Ergebnisses eine inverse Fourier- 
Transformation in den Zeitbereich durch. 

Bei diesem Ausfuhrungsbeispiel werden die in Fig. 26 gezeigten Simulationsgleichnngen der Momentenme- 
thode iediglich durch Beriicksichtigung der gegenseitigen Impedanz gelost; werden hingegen die in Fig. 41 
65 gezeigten Simulationsgleichungen verwendet, die zudem die gegenseitige Admittanz und die Wechselwirbmg 
berucksichtigen, konnen die gegenseitige Admittanz und die Wechselwirkung ebenf alls durch die zuvor erwahn- 
ten Naherungsgleichungen auBerst schnell ermittelt werden. 

Bei der in Fig. 41 dargesteHten Gleichung bezeichnet Ic,n einen in einem Metallobjekt flieBenden elektrischen 
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i • _.i j„ nhorflsche ernes dielektrischen Objekts flieBenden elektrischen Strom imd Ma 
Strom - itj^g&gS. tttSn bbiete ffieftenden magnetischen Strom. Der hochgestellte 
V^if^Sn W«rtSS2Mi hochgestellte Index "d" berieht sichauf einen Wert m dem 
SekkSLhTn Ob^dfe Evening V bezieht sich auf ein MctaQobjekt und die Erwenerung V bexieat s>ch 

^S^fdte^SmBbcispid wird nach der LSsung dieser CHeichungen der Momentenmethode das 
Auch bei ^^J^*"*^"?^ Tterechnung einer Spannung", "Berechnang ernes elek- 

%£ZSSttS*2. £2SKSS*2i Lehem 3^V— * und M-* der 

mnjjp«™"»j^ . inversen Fotmer-Transformatioii unterzogen. Die 

^^ntS'^ auf den berechneten ^berwd^ 

^SendS werden, am den entsprechenden Strom im Zeitbereich zu enmrtein, wobe, aiischlieBendun 
S^d^bengenannten Verfahrecsschritte -Berechnong einer Spaonung", •Berechnung ernes etexro- 
^gne^cnFelde^f-Smang einer zwischen Leitem anliegenden Spannung- und -Berechnung der Kopp- 

h.^rch em^n dner Antenne abgestrahhes elektrisches Feld zu untersachen(sogenannter Immtmrtatstest). 

Die^Ttr^ beLpfel^eise lie in Fig. 42 gezeigt durchgefflbn, wobe; ein TragersgnaL welches ; erne 
Fr, n ^e4 z^c^MMHz und 1GHz anfwebt. mit einer ModulationsweUe mh emer Frequenz von Ittfaand 
e^ MoSSgrad^n^ modufiert and anschlieBend die elektronische Vorrichtong mtt der moddier- 

^ oben beschriebenen Verfahren durchgefuhrt werder, wird ; tei 

SKrri&enz des Tragersignals, V c die Amplitude des TYugersignals. a* me Kreisfrequenz der Modula- 
tionsweHe and m den Modutationsgrad der ModulationsweUe bezeicnnet: 



Fonnel^V-VdCl+insiiiOmt)!^©^ 



fDahei wurde die Beziehtmg'cost — jsnrf" berucksichtigt.) , , 

T^n^ttnmrtnodewWnac^^ . „ ... fc 

Sfi^^l^e^^r^men^meAode die in dem Element m flieBenden Strome an Frequenzbereich 

tet und addiert werden, urn den in dem Element m flieBenden Strom za berechnen. . 
Slinn tefeto WellenqueOe, die eine durch Modulation ernes Trageragnab erzcugte Welle anssendet 
Somit kann planer ""j^^J^^,. ^ magnensche Strom im Zeitbereich ohne Anwendimg der 
^SfoSn eSl wSTd^ Stanch die ^em bestimm^n Elen«nt ad*re*nde 
In^uneund dasto einem bestimmten Element auftretende elekmsche und nmgnensche Feld_ ^J^O*™* 
Sr^HmSiSereich vorliegenden und wie oben dargestelh mnrttehen elel«mchen_ 
Iwme vor deren Kombination oder nach deren Kombination enmttelt werden, so daB der EinflaB der WeDen- 

^^t^^rT^^^^^^ Strom im Zeitbereicn durch Kornbuarion 
^SSSSSaSS^ "srro^ berechnet and die zwischen Leitem auliegende Spannung on Ze«bereich 
»nh^d der adchone ffir die an einem imaginar zwischen die Letter eingeffigten WWersmd autaetende 
^a^be^Sero^we" dieser Wid^and unendTich gro B gewahh ^^^^'^7°^ 
..^r^rame Iconnen vor deren Kombination ausgewertet und der in emer Empfangsschatamg etc emer 
f^SeiXauf^endTstrom und die entsprechende Spannung im Frequen^ere^h anha«ider 
^J^^S "^TavS? berechnet werden, wobei die berechneten Werte anschlieBend kombmiert werden, 
S ^rde^r^gsscSg etc de^benragun^.e^g^grund der WeflenqneDe hervorgerofene 

ana Tvon d« ^rahlnng eines von der. in Fig. 42 gezeigten Antenne erzengten elekmschen Feldes beemfhiBt 

W ^„ oben beschrieben worden ist, werden erfindungsgemaB die gegenseWge topedanz, die gegenseitige 
AdSr^ unT MMnWbd der Ausfuhrung der Momentenmethode durcb . die ™ 

werdelL wobei das Simulationsergebiiis in den Zeitbereich zurijcktransf ormiert wird- w^;^,* 
l«f Wefae kann erfinduiLgemafl im Zeitbereich der EinfluB einer WeDenqueUe airf erne elektroni^he 

SdetoS^ Vorrichtung aufgrund der RadioweUen oder des Rauschens auftreten 
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ErfmdungsgemaB kann somh auf emfache Art und Weise durch Anwendung der Momentenmethode der 
EinfhiB eincr Wellenquelle, die eine durch Modulation eines Tragersignals erzeugte Welie absendet, auf erne 
elektronische Vorrichtung im Zeitbereich simuliert werden. 

Daher kann durch die Anwendung der vorfiegenden Erfindung an einem Arbeitsplatz untersucht werden, ob 
in einer elektronischen Vorrichtung eine Fehlfunktion aufgrund eines von einer Antenne erzeugten und auf die ' 
elektronische Vorrichtung einwirkenden elektrischen Feldes auftritt oder nichL 

Patentanspruche 

1. Simulationsvorrichtung (1), welche die Momenlenmethode anwendet und eine elektronische Vorrichtung 
in Elemente aufteflt und fur vorgegebene Frequenzen eine zwischen den Hementen auftretende gegenseiti- 
ge Impedanz (Z(j) berechnet und (fie gegenseitige Impedanz sowie die den Hementen ectsprechenden 
WellenqueDen zur Simulation des in jedem Element fiieBenden Stromes gemaB der Momentenmethode 
verwendet, gekennzeichnet dorch 

Transformationsmittel (11) zum Auffinden eines Frequenzspektrums, indem eine zertliche Datenfolge einer 
Wellenquelle unterteih und darauf die Fourier-Transformation angewendet wird, und zum Vorgeben des 
gesamten Frequenzspektrums oder eines Teils des Frequenzspektrums sowie die von diesem Spektrum 
umfaBten Frequenzen als Verarbeitungsobjekt, 

Berechnungsmittel (12) Auswahlen einer Abtastfrequenz unter Berucksichtigung der von den Trans- 
formationsmitteln (11) vorgegebenen Frequenzen, Berechnen der gegenseitigen Impedanz fur die Abtast- 
frequenz, Ermittem einer Naherungsgleichung fur die gegenseitige Impedanz (Zy) anhand der berechneten 
gegenseitigen Impedanz und der Abtastfrequenz, und Berechnen der gegenseitigen Impedanz fur bestrmm- 
te von den Transformatioiismittem (1 1) vorgegebene Frequenzen anhand der Naherungsgleichung, 
Sunulationsmittel (13) zum Ermittem des in jedem der Elemente fiieBenden Stromspektrums durch Anwen- 
dung der Momentenmethode anhand der durch die Berechnungsmittel (12) berechneten gegenseitigen 
Impedfl"^ und des von den Transformationsmitteln (11) vorgegebenen Frequenzspektrums, und 
Inverse-Transfommuonsmittel (14) zum Anwenden der inversen Fourier-Transformation auf das von den 
Simulationsmitteln (13) ermittette Stromspektrum und zum Ausgeben des Ergebnisses. 

2. Simulationsvorrichtung (1), welche die Momentenmethode anwendet und eine elektronische Vorrichtung 
in Elemente aufteilt und fur vorgegebene Frequenzen eine zwischen den Hementen auftretende gegenseiti- 
ge Impedanz (Z,*j) und die gegenseitige Impedanz sowie die den Hementen entsprechenden WeQenqueUen 
zur Simulation des in jedem Element fiieBenden Stromes gemaB der Momentenmethode verwendet, ge- 
kennzeichnet durch 

Transformationsmittel (11) zum Auffinden eines Frequenzspektrums, indem eine zertliche Datenfolge einer 
Wellenquelle unterteih und darauf die Fourier-Transformation angewendet wird, und zum Vorgeben des 
gesamten Frequenzspektrums oder eines Teils des Frequenzspektrums sowie die von diesem Spektrum 
umfaBten Frequenzen als Verarb ertungsobjekt, 

Berechnungsmittel (12) Auswahlen einer Abtastfrequenz unter Berucksichtigung der von den Trans- 
formationsmittein (11) vorgegebenen Frequenzen, Berechnen der gegenseitigen Impedanz fur die Abtast- 
frequenz, Ermitteln einer Naherungsgleichung fur die gegenseitige fmprdanr (Zq) anhand der berechneten 
gegenseitigen Titi pgrfan? und der Abtastfrequenz, und Berechnen der gegenseitigen Impedanz fur bestimm- 
te von den Transfomatkmsmittem (11) vorgegebene Frequenzen anhand der Naherungsgleichung, 
Simulationsmittel (13) zum Ermittem des in jedem der Elemente fiieBenden Stromspektrums durch Anwen- 
dung der Momentenmethode anhand der durch die Berecmiungsmittel (12) berechneten gegenseitigen 
Impedanz und des von den Transformationsmitteln (11) vorgegeben Fiequenzspektrums, und, falls ge- 
wGnscht, zum gleichzeitigen Ermitteln eines Spannungsspektrums, eines elektrischen Feldspektrums und 
eines magnetischen Feldspektrums anhand des Stromspektrums und zum Auswahlen eines gewunschten 
Spektrums aus den ermhtetten Spektren als Ausgabeobjekt, und 

Inverse-Transforrnationsmhtel (14) zum Anwenden der inversen Fourier-Transformation auf das von den 
Simulationsmitteln (13) zur Ausgabe ausgewahlte Spektrum und zum Ausgeben des Ergeb nis ses 

3. Simulationsvorrichtung nach Anspruch 1 oder 2, dadurch gekennzeichnet, daB, falls das zwischen Leitern 
auftretende Spannungsspektrum zu ermitteln ist, die Simulationsmittel (13) das Spektrum der zwischen den 
Leitern aoftretenden Spannung gemaB einer B erechnungsgieichung berechnen, die aus dem von einem 
imaginar zwischen den Leitern eingesetzten Widerstand erzeugten Spannungsspektrum abgelehet ist, 
wobei der Widerstandswert unendhch gro8 gewahlt wird. 

4. Simulationsvorrichtung nach einem der vorhergehenden Anspruch e, dadurch gekennzeic hn et, daB, falls 
ein Strom- oder Spannungsspektrum an Anschlussen einer Obertragungsleitung zu ermitteln ist, die Simula- 
tionsmittel (13) die Momentenmethode ohne Berucksichtigung der Obertragungsleitung ausfuhxen, abhan- 
gig von dem Ergebnis der Momentenmethode das durch das Stromspektrum der untersuchten Elemente an 
den Anschlussen auftretende elektrische Feld ermitteln und durch Anwendung des ermittehen elektrischen 
Feldes das an den Anschlussen der Obertragungsleitung auftretende Strom- oder Spannungsspektrum 
anhand einer vorgegebenen Berechnungsgleichung berechnen. 

5. Simulationsvorrichtung nach einem der vorhergehenden Anspruche, d adur ch gekennzeichnet, 

daB, falls die Anfangsspannung eines aufgeladenen Kondensators als Wellenquelle zu untersuchen ist, die 
Transformationsmittel (11) das periocfische Signal der Anfangsspannung unterteflen und darauf die Fourier- 
Transformation anwenden, und 

daB, falls der durch eine Induktivitat flieBende Anfangsstrom als Wellenquelle zu untersuchen ist, die 
Transformationsmittel (11) das periodische Signal des Anfangsstromes unterteilen und darauf die Fourier- 
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Transformation anwenden. ■.,«,»" - ■. 

6. Simulationsvorricfatung nach einemder vorhergehenden Anspruchcdadurch gekeimzeichnet, 

daB die Berechnungsmittel (12) zusatzlich za der gegenseitigen Impedanz (Zij) auch die gegensertige 

Admittanzfy^oderdieWechsehvirta^ 

daB die Simulationsmittel (13) ihre Funktion auch unter Berujtocfatigung des magnetischen Sm>mspek- 
trums in Kombmation mit der von den Berecfanungsmitteln (12) berechneten gegenseitigen Admittanz (Y q ) 

^S^S^Sw die inverse Fourier-Transformation in Oberemsummung mit 

t%^^o^^S ] ^^^^ ^^ heDdca Anspruche^ dadurcfa gekennzeicJinet, daB die 
Berechnungsmittel (12) erne zwiscben einem Element i and einem Hement j auftretende gegensertige 
Impedanz Zij naherungsweise wie foigt berecfanen: 

Zij « e-^Rao + aif 2 + taf* + aaf 6 + 34 s + + jfbof 1 + bif + bjf 3 + bsf 5 + + ...)]■ 

wobei f die Frequenz, k die WeUenzahl, r 0 den Abstand zwischen dem Element i und dem Element j sowie j 

al2d^o^£tch Anspruch 7, dadurcb gekennzeicbnet, daB die Bere^ungsmittel (12) erne 
zwischen einem Element i und einem Element j auftretende gegensertige Admittanz naherungsweise wie 
folgtberechnen: 

y„ = c -jfaO [(ao + ai p + a2f 4 + a3f s + ^8 + . . .) + j(b^ + bif + fcf 3 + baf 5 + D4F + . . .)] 

wobei f die Frequenz, k die WeUenzahl, r 0 den Abstand zwischen dem Hement i und dem Element j sowie j 

g^im^^o^rr^Sch Anspruch 7 oder 8, dadurch gekennzeictoevdafi me^r«^gsimttel(12) 
erne zwischen einem Hement i und einem Hement j auftretende Wechselwirkung By naherungsweise wie 
foigt berecfanen: 

By = c-^Jao + atf 2 + a 2 f* + ajf 5 + a+f* + ...) + j(bof + bif 1 + W* + brf 7 + b4».+ »-)] 

wobei f die Frequenz, k cfie WeUenzahl, r 0 den Abstand zwischen dem Hement i und dem Hement j sowie j 

tT StodaSor^S" [1), wefche die Momentenmethode ^^^^^V^^^^I 
tung in Hemente aufteut und fur vorgegebene Frequenzen erne zwischen ^J^^J^^ 
ge«nseitige Impedanz (Z^) berechnet and me gegensertige Impedanz sowie die den Hementen enls^c- 
chfnden V^Uenoue^zurSimulation des in jedem Hement flieBendenStromes gemSB der Momentemne- 

SSann*^^ in drei WeUenqueDen mit umerschiedhch^ Frequenzen, 

die durch eine Tragersignalfrequenz und eine Modulaiionssignaifrequenz bestimmt werden, wobei erne 

durch Modulatioii des Tragersignals erhakene Welle als WeUenquelle verwendet wird, . 

SU^tio^S (22) zum EWtem des in jedem der Hemente ffieBenden Stromes im FrequenzbereK* 
du^we^de7^mentenmem 

Bere^ungsmhtel {73) zum Berechnen des Stromes im Zeitbereich anhand des durch die Smmlationsnnttel 
(22) ermhtelten Stromes mi Frequenzbereich. 

1 1 . SunolatKHisvorrichtung nach Anspruch 10, dadurch gckenmeichnet, Ci _ . 

daB die Simulationsmittel (22) durch Anwendnng der Momentenmethode auch den niagne^en S^un 
Frequenzbereich ermittem, wobei zusatzlich zu der gegenseitigen Impedanz (Zij) auch me gegensertige 
Admittanz (Yij) oder die Wechselwirkung (B3) oder betde berucfcnchtigt werden, und 
daB die Berechnungsmittel (23) anhand des von den Simulationsmittein (22) ermittelten magnetischen 
Stromes im Frequenzbereich auch den magnetischen Strom im Zeitbereich .berechnen. _ 
12.Simmationsverfahren, welches die Momentenmethode anwendet und erne elektromsche Vomchtong m 
Hemente aufteOt und fur vorgegebene Frequenzen eine zwiscben den Hementen auftrete^^ens^ge 
Impedanz (Zij) berechnet und gemaB der Momentenmethode den m jedem Hement ^oSttom 
anhand der gegenseitigen Impedanz sowie den den Hementen entsprecfaenden WellenqueDen sinmhert, 

luleTverf^^s^ (ST1) zum Auffinden eines Frequeiizspektrums, indern emezei^che ^tenfolge 
einer WeUenquelle unterteOt und darauf die Fourier-Transformation angewendet wird, und zum Vorgeben 
des gesamten Frequeiizspektrums oder eines Teils des Frequeiizspektrums sowie die von diesem Spektrum 
umfafiten Frequenzen als Verarbeitnngsobjeki, t^u^„ t ^A^ m 

einen zwdtenVerfahrensschritt (ST2) zum Auswahlen einer Abtastfrequenz unter Beructochtigung der m 
dem ersten Verfahrensscfaritt (ST1) vorgegebenen Frequenzen, Berechnen der gegenseitigen topedanz fin- 
die Abtastfrequenz, Ennhtein einer Nanerungsgleichung fur die gegensertige Impedanz (Zij) annand der 
berechne^egensehigen Impedanz und der Abtastfrequenz, und Berechnen der gegensemgen Impec^ 
fur bestimmte durch den ersten Verfahrensschritt (ST1) vorgegebene Frequenzen anhand der Naherungs- 

finetdrlen Verfahrensschritt (ST3) zum Ermitteln des in jedem der Hemente fheBenden^on^ekmims 
dur^h^wendung der Momentenmethode anhand der in dem zweiten Verfahrensschntt (ST2) berechneten 
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gegenseitigen Tmp^^"^ und des in dem ersten Verfahrensschritt (ST 1) vorgegeben Frequenzspektrums, 
und 

einen vierten Verfahrensschritt (ST4) zom Anwenden der inversen Fourier-Transformation auf das in dem 
dritten Verfahrensschritt (ST3) ermittelte Stromspektrnm und zum Ausgeben des Ergcbnisses. 

1 3. Simulationsverfahren, welches die Momentenmethode anwendet und eine elektroaische Vorrichtung in 
Elemente aufteflt und fur vorgegebene Frequenzen eine zwischen den Elementen auftretende gegensehige 
[mpedanz (Zij) berechnet und gemaB der Momentenmethode den in jedem Element flieBenden Strom 
anhand der gegenseitigen Impedanz sowie den den Elementen entsprechenden WeSenqueUen sim n li ert, 
gekennzeichnet durch 

einen Verfahrensschritt (ST 1) zum Auffinden ernes Frequenzspektrums, indem erne zerdiche Datenfblge 
einer WellenqueUe unterteUt und darauf die Fourier-Transformation angewendet wird, und zum Vorgeben 
des gesamten Frequenzspektrums oder eines Teils des Frequenzspektrums sowie die von diesem Spektnnn 
umfafiten Frequenzen als Verarbeitungsobjekt, 

einen zwehen Verfahrensschritt (ST2) zum Auswahlen einer Abtastfrequenz unter Berucksichtigung der in 
dem ersten Verfahrensschritt (ST1) vorgegebenen Frequenzen, Berechnen der gegenseitigen Impedanz fur 
die Abtastfrequenz, Ermitteln einer Nanerungsgleichung fur die gegensehige ImpHarrz (Zij) anhand der 
berechneten gegenseitigen Impedanz und der Abtastfrequenz, und Berechnen der gegenseitigen Impedanz 
' fur bestimmte durch den ersten Verfahrensschritt (ST1) vorgegebene Frequenzen anhand der Naherungs- 
gleichung, 

einen dritten Verfahrensschritt (ST3—ST6) zum Ermitteln des in jedem der Elemente flieBenden Strom- 
spektrums durch Anwendung der Momentenmethode anhand der in dem zwehen Verfahrensschritt (ST2) 
berechneten gegenseitigen Impedanz und des in dem ersten Verfahrensschritt (ST1) vorgegeben Frequenz- 
spektrums und, faHs gewunscht, zum gteichzeragen Ermitteln eines Spannungsspektnxms, eines elektrischen 
Feldspektrums und eines magnettschen Feldspektrums anhand des Stromspektrums und zum Auswahlen 
eines gewunschten Spektrums aus den ermittelten Spektren als Ausgabeobjekt, und 

einen vierten Verfahrensschritt zum Anwenden der inversen Fourier-Transformation auf das von den 
Sonulatransmittem (13) zur Ausgabe ausgewahlte Spektrum und zum Ausgeben des Ergebnisses. 

14. Simulationsverfahren i rarh Anspruch 12 oder 13, dadorch gekennzeichnet, 

daB in dem zwehen Verfahrensschritt (ST2) zusatzlich zu der gegenseitigen Imp erfa n z (Zij) auch die 
gegensehige Admittanz (Yy) oder die Wechserwirkung (Bij) oder beide durch daBelbe Verfahren berechnet 
werden, ' 

daB der dritte Verfahrensschritt (ST3) auch unter Berucksichtigung des magnetischen Stromspektrums in 
{Combination mrt der von den Berectaungsmitteln (12) berechneten gegenseitigen Admittanz (Yq) und 
Wechseftvirkung (Bij) durchgefuhrt wird, and 

dafi in dem vierten Verfahrensschritt (ST4) die inverse Fourier-Transformation in Obereinstimmimg mit 
dem dritten Verfahrensschritt (ST3) durdigefuhrt wird 

15. Siniulationsverfahren, welches die Momentenmethode anwendet und eine elektroaische Vorrichtung in 
Elemente aafteflt und fur vorgegebene Frequenzen eine zwischen den Elementen auftretende gegensehige 
Impedanz (Zij) berechnet und den in jedem Element flieBenden Strom gemaB der Momentenmethode 
anhnnri der gegenseitigen Impedanz sowie den den Elementen entsprechenden WeOenqueHen simuliert, 
gekennzeichnet durch 

einen ersten Verfahrensschritt zum Aufteflen einer WeUenquelle in drei WellenqueBen mit unterschiedli- 
chen Frequenzen, die durch eine Tragersignalfrequenz und eine Modulatkwrssigna Ifrequenz bestimmt sind, 
wobei eine durch Modulation des Tragersignals erhaltene Welle als die WeOenqueUe verwendet wird, 
einen zwehen Verfahrensschritt zum Errnittetn des in jedem Element flieBenden Stromes im Frequenzbe- 
reich durch Anwenden der Momentenmethode auf die in dem ersten Verfahrensschritt aufgeteOten Wellen- 
queUen,und 

einen dritten Verfahrensschritt zum Berechnen des Stromes im Zeitbereich anhand des in dem zwerten 
Verfahrensschritt ermittelten Stromes im Frequenzbereich. 

16. Simulationsverfahren nach Anspruch 15, dadurch gekennzeichnet, 

daB in dem zwehen Verfahrensschritt durch Anwendung der Momentenmethode auch der magnetischen 
Strom im Frequenzbereich ermittelt wird, wobei zusatzlich zu der gegenseitigen TmpHan7 (Zq) auch <fie 
gegensehige Admittanz (Yij) oder die Wechserwirkung (Bq) oder beide berucksichtigt werden, und 
daB in dem dritten Verfahrensschritt auch der magnetische Strom im Zeitbereich anhand des in dem 
zwehen Verfahrensschritt ermittelten magnetischen Stromes im Frequenzbereich berechnet wird. 
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Fig.44A< 



a m i+jbmi= I mi exp(j <t> i) 

tan d> i=b m i/a m i 

a m 2+jbm2= I m 2 exp(j d>z) 

tan ^2 = bm2/a m 2 

am3+jbm3= I m3 eXp(j & 3 ) 

tan <£3 = b„, 3 /a m 3 



I m tsin(ojc t+<£ i) 
Fi g.AAB\ Im2sin[(<y«-Gj..)t+^2] 

I m3sin[(cuc+ajm)t+a!)3] 
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Fig. 45 
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